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ABSTRACT 

We  have  developed  a  numerical  model  for  the  temporal  evolution  of  particle  and  photon  spectra  resulting 
from  nonthermal  processes  at  the  shock  fronts  formed  in  merging  clusters  of  galaxies.  Fermi  acceleration  is 
approximated  by  injecting  power-law  distributions  of  particles  during  a  merger  event,  subject  to  constraints 
on  maximum  particle  energies.  We  consider  synchrotron,  bremsstrahlung,  Compton,  and  Coulomb 
processes  for  the  electrons,  nuclear,  photomeson,  and  Coulomb  processes  for  the  protons,  and  knock-on 
electron  production  during  the  merging  process.  The  broadband  radio  through  7-ray  emission  radiated  by 
nonthermal  protons  and  primary  and  secondary  electrons  is  calculated  both  during  and  after  the  merger 
event.  Using  ROSA  T  observations  to  establish  typical  parameters  for  the  matter  density  profile  of  clusters  of 
galaxies,  we  find  that  typical  merger  shocks  are  weak  and  accelerate  particles  with  relatively  soft  spectra.  We 
consider  the  prospects  for  detecting  nonthermal  radio  and  7-ray  emission  from  clusters  of  galaxies  and  impli¬ 
cations  for  the  origin  of  ultra-high-energy  cosmic  rays  and  the  diffuse  7-ray  background.  Our  results  suggest 
that  only  a  few  of  the  isotropically  distributed  unidentified  EGRET  sources  are  due  to  shocks  formed  in 
cluster  mergers  and  that  only  a  minor  contribution  to  the  diffuse  extragalactic  7-ray  background  can 
originate  from  cluster  merger  shocks.  Cluster  merger  shocks  can  accelerate  protons  to  <  1019  eV  for  the 
standard  parameters  considered  here.  We  predict  that  GLAST  will  detect  several  cluster  mergers,  and 
depending  on  the  mean  magnetic  fields  in  the  intracluster  medium,  the  Low  Frequency  Array  could  detect 
anywhere  from  several  to  several  hundred. 

Subject  headings:  acceleration  of  particles  —  galaxies:  clusters:  general  —  methods:  numerical  — 
shock  waves  —  X-rays:  galaxies:  clusters 


1.  INTRODUCTION 

According  to  the  hierarchical  merging  scenario,  cold  dark 
matter  halos  evolve  to  form  larger  structures  by  merging 
with  adjacent  dark  matter  halos.  Within  dark  matter  halos, 
baryonic  matter  condenses  to  form  clusters  of  galaxies.  A 
cluster  merger  event  results  from  the  interaction  of  galaxy 
clusters  during  the  merger  of  cold  dark  matter  halos.  The 
gravitational  potential  energy  available  in  a  cluster  merger 
event  involving  halos  with  masses  ~1015  M0  is  ~1063-1064 
ergs.  As  two  clusters  of  galaxies  merge,  the  infall  velocities 
can  exceed  the  sound  speed  of  the  intracluster  medium 
(ICM).  As  a  result,  a  shock  front  will  form  at  the  interaction 
boundary  between  the  clusters.  First-order  Fermi  accelera¬ 
tion  at  the  shock  front  produces  a  population  of  non¬ 
thermal,  relativistic  particles.  Relativistic  electrons  are 
detected  from  their  synchrotron  radio  emission  or  from 
7-rays  due  to  Compton-scattered  cosmic  microwave  radia¬ 
tion.  Nonthermal  protons  are  detected  through  7-rays  emit¬ 
ted  from  secondaries  formed  in  nuclear  production 
processes,  including  the  7r°-decay  signature  at  70  MeV. 

Rich  clusters  contain  a  hot  and  tenuous  ionized  ICM, 
with  observed  temperatures  Tx  ~  5-10  keV,  sound  speeds 
~1000  km  s  1 ,  and  thermal  bremsstrahlung  luminosities 
Lx  ~  1045  ergs  s  1  between  2  and  10  keV.  In  addition  to  the 
luminous  thermal  component  present  in  these  clusters,  there 
is  a  growing  body  of  evidence  supporting  the  presence  of 
nonthermal  distributions  of  particles  in  cluster  mergers 
(Eilek  1999;  Feretti  et  al.  2000).  Deep  radio  observations  of 
clusters  of  galaxies  indicate  the  presence  of  extended  diffuse 
emission  not  easily  associated  with  an  optical  counterpart. 
These  diffuse  radio  features  are  commonly  classified  as 
either  radio  halos  or  radio  relics  (which  are  also  called 


periphery  halos).  Radio  halos  mimic  the  observed  X-ray 
profiles  and  are  characterized  by  their  central  location  in  the 
cluster  and  by  a  highly  disorganized  magnetic  field  and  are 
generally  thought  to  be  a  consequence  of  a  merger  event 
(Feretti  et  al.  2000).  Radio  relics,  usually  found  on  the 
periphery  of  the  cluster,  are  characterized  by  highly  organ¬ 
ized  magnetic  fields  and  often  display  filamentary  struc¬ 
tures.  We  focus  here  on  radio  relics,  which  are  thought  to 
result  from  synchrotron  emission  emitted  by  electrons 
directly  accelerated  at  shock  fronts.  By  contrast,  radio  halos 
might  be  due  to  reacceleration  of  relic  electrons  by  magnetic 
turbulence  or  enhanced  magnetic  turbulence  arising  from 
the  cluster  merger  or  motions  of  the  galaxies  (Brunetti  et  al. 
2001;  Ohno,  Takizawa,  &  Shibata  2002). 

The  first  radio  halo  in  a  cluster  of  galaxies  was  detected 
from  the  Coma  Cluster  (Large,  Mathewson,  &  Haslam 
1959).  Many  halos  and  relics  have  since  been  found  in  a 
number  of  other  clusters,  including  A754  (Kassim  et  al. 
2001),  A2256  (Bridle  &  Fomalont  1976),  and  others  (see 
Govoni  et  al.  2001;  Slee  et  al.  2001  for  recent  observations). 
Typical  radio  powers  from  the  radio  halos  and  relics  are  at  a 
level  of  1040  1  042  ergs  s  1  (Giovannini  &  Feretti  2000). 
These  clusters  show  a  favorable  correlation  (Kassim  et  al. 
2001;  Berrington,  Lugger,  &  Cohn  2002)  between  the  exis¬ 
tence  of  diffuse  radio  features  and  recent  or  ongoing  merger 
activity.  Detection  of  these  features  is  largely  due  to  sensitiv¬ 
ity  limitations;  the  number  of  detected  diffuse  radio  features 
has  increased  with  the  improvement  in  sensitivity  of  radio 
telescopes. 

The  synchrotron  evidence  for  nonthermal  electrons  in 
clusters  of  galaxies  indicates  that  an  acceleration  mecha¬ 
nism  must  exist  in  the  environment  of  the  host  cluster  with 
sufficient  power  to  produce  the  observed  emission.  Several 
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mechanisms  have  been  proposed  to  explain  the  correlation 
of  the  radio  halo  and  relic  features  and  the  recent  or  ongoing 
merger  activity.  These  mechanisms  often  require  the  pres¬ 
ence  of  shock  waves  or  magnetic  turbulence  to  accelerate 
particles  via  the  Fermi  acceleration  process  (Schlickeiser, 
Sievers,  &  Thiemann  1987;  Tribble  1993;  Kang,  Rachen,  & 
Biermann  1997;  Ensslin  et  al.  1998;  Blasi  2000;  Miniati  et  al. 
2001).  Other  theories  include  adiabatic  compression  of 
fossil  radio  plasma  by  a  cluster  merger  shock  wave  (Ensslin 
&  Gopal-Krishna  2001). 

Optical  surveys  show  that  approximately  30%-40%  of 
clusters  of  galaxies  display  evidence  for  the  presence  of  sub¬ 
structure  (Forman  et  al.  1981;  Geller  &  Beers  1982  and 
others).  This  internal  structure  is  often  interpreted  as  a  sub¬ 
set  of  galaxies  merging  with  a  larger  cluster  of  galaxies. 
These  internal  structures  indicate  velocity  differences  near 
or  greater  than  the  expected  sound  speed  of  the  ICM. 
Observed  velocities  typically  range  from  ~  1 000  to  3000  km 
s_I,  which  are  consistent  with  values  expected  from 
parabolic  orbits  (Oegerle  &  Hill  1994). 

Emission  from  nonthermal  particles  will  also  appear  at 
EUV  and  hard  X-ray  (HXR)  energies  as  a  power-law  excess. 
While  the  EUV  emission  seen  in  the  Coma  Cluster  (Lieu 
et  al.  1999b)  and  possibly  also  A2199  and  A1795  (Lieu, 
Bonamente,  &  Mittaz  1999a)  may  have  a  cool  thermal  ori¬ 
gin,  it  is  unlikely  because  of  the  extreme  mass  requirements 
of  cool  gas.  It  is  more  likely  that  the  EUV  emission  has  a 
nonthermal  origin  (Hwang  1997;  Ensslin  &  Biermann  1998; 
Mittaz,  Lieu,  &  Lockman  1998;  Blasi  &  Colafrancesco 
1999;  Atoyan  &  Volk  2000).  Excess  HXR  emission  has  been 
reported  from  the  clusters  A1656  (Fusco-Femiano  et  al. 
1999;  Rephaeli,  Gruber,  &  Blanco  1999),  A2256 
(Fusco-Femiano  et  al.  2000),  A3667  (Fusco-Femiano  et  al. 
2001),  and  possibly  A2199  (Kaastra  et  al.  1999).  Estimated 
luminosity  of  the  HXR  emission  is  «1043  ergs  s  1 .  Thermal 
origins  of  this  HXR  excess  emission  would  require  unrealis¬ 
tic  temperatures  of  greater  than  40  keV,  and  so  it  is  thought 
to  be  caused  by  the  presence  of  relativistic  nonthermal 
electrons  (Fusco-Femiano  et  al.  1999). 

Numerical  models  of  merging  clusters  of  galaxies  (Ricker 
1998;  Takizawa  2000;  Ricker  &  Sarazin  2001;  Miniati  et  al. 
2001  and  others)  have  treated  the  development  of  shocks  as 
a  result  of  the  cluster  merging  process.  Given  the  presence 
of  thermal  ionized  particles  in  the  vicinity  of  these  shocks 
and  a  cluster  magnetic  field,  electrons  and  ions  will  be  accel¬ 
erated  via  the  first-order  Fermi  process  (Bell  1978; 
Blandford  &  Ostriker  1978).  Numerous  attempts  have  been 
made  to  model  the  emissions  from  nonthermal  particles 
(Colafrancesco  &  Blasi  1998;  Fujita  &  Sarazin  2001; 
Petrosian  2001;  Miniati  et  al.  2001)  produced  by  these  clus¬ 
ter  merger  shocks,  but  this  study  differs  from  previous 
attempts  in  that  we  accurately  model  the  diffusion  of  par¬ 
ticles  in  energy  space  due  to  Coulomb  interactions  and 
allow  for  a  variable  injection  rate  that  depends  on  the  envi¬ 
ronment  local  to  the  shock  front.  Recent  studies  by  Liang, 
Dogiel,  &  Birkinshaw  (2002)  have  addressed  the  balance  of 
Coulomb  losses  and  diffusion  in  energy  space  but  do  not 
address  the  injection  and  diffusion  of  particles  in  energy 
space  along  with  a  variable  source  function  of  nonthermal 
particles.  In  our  treatment,  we  follow  particle  energies  up  to 
1019  eV  and  calculate  the  bremsstrahlung,  Compton,  syn¬ 
chrotron,  and  7ro  7-rays  from  p-p  collisions.  Because  the 
shock  front  lifetime  is  a  significant  fraction  of  the  age  of  the 
universe,  we  also  follow  the  changing  environment  due,  for 


example,  to  the  changing  cosmic  microwave  background 
(CMB)  energy  density.  In  addition,  we  accurately  model 
nonthermal  electrons  and  protons  up  to  ~1021  eV,  although 
we  find  that  limitations  on  particle  acceleration  make  it  diffi¬ 
cult  to  produce  protons  above  ~  1 019  eV  in  cluster  merger 
shocks.  From  the  particle  distributions,  we  calculate  non- 
thermal  photon  spectra  for  energies  up  to  ~107  GeV.  We 
also  include  the  effects  of  secondary  production  on  the 
nonthermal  photon  spectra. 

The  physical  processes  and  the  temporal  evolution  of  the 
particles  are  described  in  §  2.  The  results  of  the  simulations 
are  presented  in  §  3.  A  comparison  of  the  photon  spectra 
with  observed  clusters  and  the  potential  of  detecting  these 
shocks  with  space-based  satellite  observatories  is  discussed 
in  §4. 


2.  MODELS 

We  have  developed  a  code  to  calculate  the  time-depend¬ 
ent  particle  distribution  functions  evolving  through  radia¬ 
tive  losses  for  electrons  and  protons  accelerated  by  the  first- 
order  Fermi  process  at  the  cluster  merger  shock.  The  code 
was  originally  adapted  from  a  supernova  remnant  code 
(Sturner  et  al.  1997)  and  applied  to  the  specific  case  of  the 
cluster  merger  scenario.  We  break  up  the  discussion  of  the 
problem  into  four  subsections:  cluster  merger  dynamics, 
nonthermal  particle  production,  temporal  evolution,  and 
the  production  of  the  photons. 


2. 1 .  Cluster  Merger  Dynamics 
2.1.1.  Evolution  of  the  Merger  Event 

Two  clusters  are  assumed  to  be  the  dominant  gravitating 
masses  in  the  local  region,  and  variations  of  the  gravita¬ 
tional  field  due  to  surrounding  masses  are  assumed  to  be 
negligible.  A  cluster  accretes  another  cluster  that  falls  from 
a  distance  d.  This  distance  can  be  approximated  by  the  max¬ 
imum  separation  of  two  point  masses  that  move  apart 
beginning  at  the  moment  of  the  big  bang,  follow  an  elliptical 
orbit,  and  merge  at  time  tm.  Given  this  assumption,  a  cluster 
of  mass  M\  =  1015Mi5  AT.  accretes  a  subcluster  of  mass 
Mi  =  1015/«i5  Mq  at  time  tm,  as  approximated  by  Kepler’s 
third  law.  Consequently,  a  cluster  that  merges  at  a  time  t,„ 
fell  from  a  maximum  separation  approximated  by 


2G(M\  +  M2)tU  1/3 


A  M\+  M2\ 1/3  (  tm  \2'\ , 

' 45  \  1015  M&)  (lO'OyrJ  MpC  (1) 


(Ricker  &  Sarazin  2001).  This  equation  allows  us  to  calcu¬ 
late  the  magnitude  of  the  energy  pool  available  for  particle 
acceleration.  The  total  orbital  energy  §  of  the  merging 
system,  assuming  a  zero  impact  parameter,  is  given  by 


M\  M2 


21/3(M1  +  M2)1/3 

Mentis 


ttG\ 

W 


2/3 


-1.9  x  10 


64 


(Mis  +ffli5)1/3  V1010  yr 


-2/3 


ergs  .  (2) 


We  use  the  terminology  that  orbits  are  bound  if  the  total 
energy  $  <  0;  otherwise,  they  are  unbound.  It  is  customary 
to  use  the  Keplerian  analogy  of  elliptical  to  represent  orbits 
where  &  <  0  and  parabolic  for  orbits  where  $  =  0. 
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To  approximate  the  merger  velocity  of  the  system,  we 
relax  the  assumption  that  the  accreting  cluster  is  a  point 
mass.  Instead,  we  assume  the  ICM  number  density  of  the 
accreting  cluster  follows  the  spherically  symmetric, 
isothermal  beta  model  given  by 


»icmO)  =  no 


1  + 


(3) 


The  quantity  n0  is  the  central  number  density  of  the  ICM, 
and  rc  is  the  core  radius,  which  characterizes  the  spatial  scale 
over  which  the  density  changes  from  a  roughly  constant 
value  at  r<rc  to  a  power-law  density  dependence  at  r>rc. 
The  gravitational  dynamics  are  dominated  by  the  dark  mat¬ 
ter  distribution  of  the  more  massive  cluster.  Assuming  that 
the  dark  matter  distribution  follows  a  profile  similar  to  that 
of  the  ICM  and  that  the  merging  cluster  is  a  point  mass,  we 
calculate  the  cluster  merger  velocity  and  radial  separation 
by  solving  the  differential  equation 

d2rn  d$[Mi{rn),rn] 

=  — •  (4) 


The  radius  r]2  is  the  separation  of  the  center  of  masses  of  the 
two  clusters,  is  the  mass  interior  to  the  subcluster  at 

radius  r]2,  M2  is  the  mass  of  the  merging  subcluster,  and 
mr  =  (M\M2)/{M\  +  M2)  is  the  reduced  mass.  The 
gravitational  potential  <b(Mi ,  r\2)  is  defined  by 


GM{r ) 
r 


(5) 


The  mass  interior  is  defined  as 


M(r)  =  4<m  [  dxx2 hicmW  ,  (6) 

Jo 

where  is  a  normalization  constant  chosen  so  that  M\  (ri2) 
is  equivalent  to  the  total  mass  of  the  dominant  cluster  at  its 
maximum  radius  R\ .  Obviously,  the  potential  satisfies  the 
requirement  <b(M,  r  — >  oo)  =  0.  Initial  conditions  (t  =  0) 
are  found  by  a  conservation  of  energy  argument  where  the 
initial  radius  and  infall  velocity  are  given  by 

r\i(t  =  0)  =  R\  +  R2  , 


vo(t  =  0)  =  <  —  [S  —  r\2)\  ,  (7) 

y  nir 

where  R2  is  the  maximum  radius  of  cluster  M2. 

In  order  to  test  the  accuracy  of  this  approach  to  the 
dynamics  of  cluster  mergers,  we  compare  the  semianalytic 
values  of  r\2  and  vo  with  the  results  of  an  A'-body  simulation 
of  a  two-cluster  merger  model  (Berrington  2000).  For  this 
comparison,  we  assume  that  the  two  clusters  have  equal 
masses.  This  is  the  worst  case,  because  when  one  cluster  is 
smaller  and  less  massive  than  the  other,  it  is  better  described 
as  a  point  source.  In  the  A'-body  simulation,  each  cluster  is 
approximated  by  N  =  80, 000  equal-mass  particles,  with 
50%  of  the  particles  assigned  to  galaxies  in  each  cluster  and 
the  remaining  50%  assigned  to  a  nonluminous  dark  matter 
that  follows  a  King-model  distribution  with  Wo  =  6.25 
(King  1966).  The  kinetic  energy  of  the  particles  in  each  clus¬ 
ter  is  assumed  to  be  in  virial  equilibrium  with  respect  to  the 
total  gravitational  potential  energy  prior  to  the  interaction. 


D 


Fig.  1. — Comparison  between  the  semianalytic  model  used  in 
subsequent  calculations  and  results  of  an  N-body  simulation.  The  initial 
configuration  of  an  IV-body  simulation  of  a  cluster  merger  between  two 
equal-mass  clusters  is  shown  in  the  top  panel.  The  relative  differences  Al¬ 
and  An  between  the  semianalytic  and  numerical  calculations  of  the  center- 
of-mass  distance  and  the  infall  velocity  of  the  two  clusters  are  shown  in  the 
bottom  panel. 


Each  cluster  consists  of  50  galaxies  with  masses  that  are 
selected  from  a  Schechter  luminosity  function  with  a  power- 
law  index  of  —1.25.  The  top  panel  of  Figure  1  shows  the 
initial  configuration  of  the  merger  event,  and  the  bottom 
panel  gives  the  relative  difference  between  the  centers  of 
mass  calculated  in  the  A'-body  simulation  and  semi- 
analytically. 

The  calculation  ends  at  the  collision  time  rcou,  which  is 
defined  as  the  time  when  the  centers  of  mass  of  the  two  clus¬ 
ters  coincide.  As  can  be  seen,  the  semianalytic  value  of  Po  is 
accurate  to  within  10%  at  times  preceding  fcou .  After  fcon, 
the  forward  shock  velocity,  v\,  can  no  longer  be  approxi¬ 
mated  by  equation  (7)  because  the  shock  decouples  from  the 
gravitational  infall  of  the  merging  cluster  and  is  assumed  to 
be  constant  as  the  shock  propagates  outward  from  the 
cluster  center. 

In  this  paper,  we  model  the  cluster  merger  shocks  rather 
than  the  accretion  shocks.  An  infalling  cluster  is  assumed  to 
be  a  virialized,  quasi-spherical  distribution  of  galaxies, 
ICM,  and  dark  matter.  A  cluster  merger  shock  forms  when 
the  cluster  merger  velocity  exceeds  the  sound  speed  of  the 
ICM.  In  contrast,  an  accretion  shock  is  the  result  of  a 
spherically  symmetric  infall  of  matter  from  the  surrounding 
volume  onto  a  cluster.  In  the  context  of  the  previous  work 
by  Fujita  &  Sarazin  (2001),  our  definition  of  merger  events 
includes  not  only  their  definition  of  a  cluster  merger  event 
(M2  >  O.6M1 )  but  also  their  definition  of  semimerger  events 
(0.6M\  >  M2  >  0.1  Mi).  In  their  calculations,  accretion 
events  can  account  for  only  «10%  of  the  nonthermal  pho¬ 
ton  luminosity.  Although  each  cluster  of  galaxies  is  at  the 
center  of  an  accretion  flow,  we  do  not  consider  nonthermal 
particle  production  from  accretion  shocks  because  their 
nonthermal  photon  production  is  negligible  in  comparison 
to  cluster  merger  events.  Note  that  supersonic  cluster 
merger  shocks  are  less  likely  to  form  at  high  redshift  for  two 
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clusters  of  galaxies  with  given  masses.  This  is  because  the 
relative  velocities  of  the  merging  clusters  will  be  associated 
with  orbits  with  smaller  separations  at  earlier  times  and  will 
therefore  have  smaller  relative  velocities.  In  reality,  the 
mean  cluster  masses  are  smaller  at  early  times.  These  lower 
mass  clusters  have  a  lower  virial  temperature.  If  this  is  con¬ 
sidered  in  the  calculation  of  the  strengths  of  the  shocks  seen 
in  cluster  merger  events,  it  is  possible  that  stronger  shocks 
may  be  seen  at  higher  redshifts. 


where  x=  f 2(z)  -  1,  and  H(z)  =  £70(  1  +  z)3 /E2(z).  These 
relations  are  accurate  to  1%  in  the  range  fi(z)  =  0. 1-1 .0. 

A  number  of  radiation  fields  provide  seed  photons  for  the 
Compton  scattering  by  the  nonthermal  electrons  present  in 
the  cluster.  For  the  X-ray  radiation  field,  we  calculate  the 
X-ray  luminosity  of  a  cluster  of  mass  M  by  the  observed 
relation  (Arnaud  &  Evrard  1999) 


Lx(Tx)  =  2.88/C2 


x  1044 


'kTx(z) 
6  keV 


2.88 

ergs  s 


(12) 


2. 1 .2.  The  Cluster  Environment 

Shock  formation  and  nonthermal  particle  production 
depend  on  the  properties  of  ICM,  which  are  largely  deter¬ 
mined  by  the  mass  of  the  cluster  and  cosmological  epoch. 
Below  we  describe  the  equations  used  to  describe  the 
properties  of  the  cluster  environment. 

X-ray  observations  of  galaxy  clusters  show  that  these 
clusters  are  permeated  by  a  hot  ICM  of  temperatures  in  the 
range  of  5-12  keV.  At  these  temperatures  the  gas  is  well 
described  by  an  ideal  gas.  The  sound  speed  of  an  ideal  gas  in 
terms  of  the  gas  temperature  kTx  is  simply 

MTx)  =  /H  =  12“'r',! (tUv)  /  tm  s_l  ■  (8) 

where  F  represents  the  ratio  of  specific  heats,  which  we  take 
to  be  equal  to  5/3,  /u  is  the  mean  atomic  mass,  and  mp  is  the 
proton  mass. 

The  temperature  of  the  ICM  depends  on  both  the  mass  of 
the  cluster  and  redshift  z.  In  order  to  calculate  the  cluster 
X-ray  luminosity,  which  gives  the  energy  of  X  rays  that  can 
be  Compton  scattered,  we  first  calculate  the  expected  tem¬ 
perature  of  the  ICM  gas.  We  use  the  cluster  M-T  relation¬ 
ship  from  Bryan  &  Norman  (1998),  given  by 


where  kTx  is  defined  by  equation  (9).  Assuming  the  X-ray 
luminosity  is  uniformly  distributed  within  a  spherical  cluster 
of  radius  R  and  an  escape  time  of  ~R/c,  the  mean  energy 
density  of  the  X-ray  photons  is 


Ux(Tx) 


1.5  h1  x  10 


-10 


'kTx(z) 
6  keV 


-Ti — Y 

1  Mpc ) 
MeV  cm 


-3 


(13) 


The  thermal  X-ray  bremsstrahlung  is  approximated  by  a 
blackbody  of  temperature  Tx  in  our  calculations  of 
Compton  scattering.  This  is  not  a  crucial  approximation, 
because  the  X-ray  energy  density  is  ~2-3  orders  of  magni¬ 
tude  less  than  the  CMB  energy  density,  and  the  importance 
of  Compton-scattered  X-ray  photons  is  furthermore 
reduced  because  of  Klein-Nishina  effects. 

A  second  component  that  contributes  to  the  mean 
photon  energy  density  of  the  ICM  is  the  stellar  photon  field. 
We  assume  that  the  galaxy  luminosity  function  is  well 
approximated  by  the  Schechter  luminosity  distribution: 

i(£)  <>4> 


kTx(z)  =  1.39 '  [Ach2E2(z))1/3  keV  ,  (9) 

where  M\  is  the  mass  of  the  cluster,  h  is  the  parameter 
defined  by  H  =  100  /z  km  s-1  Mpc-1,  the  function  E2  (z)  = 
O0(1  +  zy  +  LlR(l  +  z)2  +  UA,  and  fTx  is  a  normalization 
constant  taken  to  be  0.8.  We  have  made  use  of  the  following 
standard  definitions: 

«o=^§?,  n*  =  w*r2,  (10) 

where  H0  is  the  current  measured  Hubble  constant,  pc  is  the 
critical  density  required  to  close  the  universe,  A  is  the  cos¬ 
mological  constant,  and  R  is  the  current  curvature  radius  of 
the  universe.  We  force  the  constraint  Q0  +  QR  +  UA  =  1  as 
mandated  by  an  inflationary  cosmology.  For  all  models  in 
this  paper,  we  set  flR  =  0  for  a  flat  universe.  The  function  Ac 
is  the  critical  density  factor  and  marks  the  point  at  which  a 
region  of  overdensity  makes  the  transition  from  linear  den¬ 
sity  growth  behavior  to  a  nonlinear  density  growth  regime. 
This  marks  the  birth  of  a  cluster  or  the  beginning  of  the 
initial  collapse  of  a  density  perturbation  out  of  the  Hubble 
flow.  Bryan  &  Norman  (1998)  calculate  it  to  be 

r  18tt2  ,  Ho  =  1  , 

Ac  =  l  18tt2  +  82x  -  39x2  ,  H*  =  0  ,  (11) 

[  18tt2  +  60.v-32x2  ,  HA  =  0  , 


where  «gai  ( L)  dL  is  the  differential  luminosity  function  and 
is  defined  as  the  number  of  galaxies  within  the  luminosity 
range  L  and  L  +  dL.  The  parameter  L*  is  the  “  characteris¬ 
tic  luminosity”  of  a  galaxy  that  marks  the  transition 
between  the  faint-end  power  law  of  slope  - a  and  the  high- 
end  exponential  cutoff.  The  constant  £*  normalizes  the 
Schechter  function  to  /Vgai  when  integrated  over  all  luminos¬ 
ities.  We  can  use  equation  (14)  to  calculate  the  total  stellar 
luminosity  of 

/»oo 

Lstar  =  /  dL  Lnga\(L)  =  £*r(2  -  o)L*  ,  (15) 

Jo 

where  T(2  -  a)  is  the  gamma  function.  Typical  values  for 
£*,  L*,  and  a  are  given,  e.g.,  in  Paolillo  et  al.  (2001),  and  our 
adopted  values  are  a  =  1.1,  a  value  of  L*  ^  7  x  1043  ergs 
s-1  corresponding  to  an  absolute  magnitude  M*  =  —20.5 
in  Bj  (Colless  1989),  and  £*  =  100. 

Using  the  same  approach  as  for  the  X-ray  photon  field, 
the  estimated  mean  stellar  radiation  field  becomes 

Ustar  M  1.4/;2  x  10~9  ( T-^—]  MeV  cnU3  .  (16) 

\1  Mpc ) 

The  spectral  distribution  of  elliptical  galaxies  is  well 
approximated  by  a  K3  III  star  (Pierce  &  Berrington  2003). 
Since  higher  density  regions  in  rich  galaxy  clusters  are 
dominated  by  elliptical  galaxies,  we  assume  that  the  stellar 
radiation  energy  distribution  can  be  approximated  by  a 
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blackbody  spectrum  resembling  that  of  a  K3  III  star  at 
temperature  7’star  ~  4000  K. 

The  energy  density  of  the  CMB  is 

C/cmb  =  2.5  x  10~7(1  +z)4  MeV  cm'3  .  (17) 

The  CMB  is  the  dominant  contributor  to  total  cluster 
photon  energy  density  by  approximately  3  orders  of  magni¬ 
tude.  The  total  cluster  radiation  field  is  simply  the  sum  of 
equations  (13),  (16),  and  (17). 

2.1.3.  Cluster  Redshift 

The  cluster  environment  described  in  §  2.1.2  depends  on 
the  redshift  of  the  cluster.  The  redshift  evolution  of  the  clus¬ 
ter  depends  on  the  chosen  cosmological  model.  We  have 
chosen  to  uniquely  specify  the  cosmological  model  by  the 
following  parameters:  the  critical  mass  fraction  f l0,  the 
curvature  ClR,  the  dark  energy  fiA,  and  the  Hubble  constant 
H0.  For  this  paper  we  consider  the  model  (fl0,  Hr,  Ha)  = 
(0.3,  0.0,  0.7). 

The  observed  redshift  of  a  cluster  is  determined  from 
the  time  since  the  big  bang  by  the  cosmological  equation 
(Peebles  1993) 

jt  =  -H0(l+z)E(z),  (18) 

where  E2(z)  is  the  dimensionless  parameter  defined  follow¬ 
ing  equation  (9).  Typical  propagation  times  for  a  shock 
front  to  traverse  a  cluster  are  ~109  yr.  This  is  a  significant 
time  span  in  that  the  evolution  of  the  cluster  environment 
due  to  cosmic  expansion  must  be  considered.  Cluster  merger 
events  can  be  initiated  at  any  time  throughout  age  of  the 
universe,  and  evolving  the  cluster  environment  is  necessary 
to  accurately  calculate  the  nonthermal  photon  and  particle 
spectra. 


2.2.  Nonthermal  Particle  Production 

Henceforth,  we  denote  particle  kinetic  energies  by 
Kep  =  mej>( lej,  —  1  )c2,  total  energies  by  Eep,  and  dimen¬ 
sionless  total  energies  by  e,-,  with  the  subscripts  e  and  p  refer¬ 
ring  to  electrons  and  protons,  respectively.  First-order 
Fermi  acceleration  at  a  cluster  merger  shock  is  approxi¬ 
mated  by  injecting  power-law  momentum  spectra  for  the 
electron  and  protons.  The  total  particle  injection  function  in 
terms  of  kinetic  energy  is  given  by 

Qe,P(Ke^  t )  =  Q% [Kej,(Kej,  +  2 mej,c2)\ ~W+l]/2 


£max(0j  ’ 

where  s(t)  is  the  injection  index  and  Amax  is  the  maximum 
particle  energy  determined  by  three  conditions:  the  avail¬ 
able  time  to  accelerate  to  a  given  energy  since  the  beginning 
of  the  merger  event,  the  requirement  that  the  particle 
Larmor  radius  is  smaller  than  the  size  scale  of  the  system, 
and  the  condition  that  the  energy-gain  rate  through  first- 
order  Fermi  acceleration  is  larger  than  the  energy-loss  rate 
due  to  synchrotron  and  Compton  processes.  The  duration 
of  particle  injection  is  determined  by  the  dynamics  of  the 
system,  and  the  shock  infall  velocity  is  determined  by  equa¬ 
tion  (4).  The  shocks  terminate  at  the  time  taCc,i  for  the  for¬ 
ward  shock  and  tacc2  for  the  reverse  shock  when  the  shock 
propagates  to  the  outer  part  of  the  cluster.  The  time  tacc 


x  (Kej>  +  mepc2)  exp 


represents  the  greater  value  of  two  acceleration  times  4cc,i 
and  tacc, 2-  The  constant  Q°e,p  normalizes  the  injected  particle 
spectrum  over  the  volume  swept  out  by  the  shock  front  and 
is  determined  by 


dKe  p  Kep  Qr  .P  ( Ec.p ,  i  ) 


^T^7?He"VH(03 («ICm(0)  , 

(20) 


where  A  is  the  area  of  the  shock  front,  rj^e  =  1 .2  is  an 
enhancement  factor  due  to  the  presence  of  ions  heavier  than 
hydrogen,  and  («icm(0)  is  the  number  density  of  the  gas 
averaged  over  the  area  of  the  shock  front  at  time  t.  We 
assume  a  constant  efficiency  factor  r/ap  =  5%  for  both  elec¬ 
trons  and  protons,  although  rjej,  would  depend  on  Mach 
number  through  the  uncertain  wave  generation  processes  at 
the  shock  front.1  The  particle  injection  index,  s(t),  is  deter¬ 
mined  by  the  shock  properties  described  in  §  2.2.1.  Typical 
values  for  the  particle  injection  index  are  s  «  2.5. 

The  temporally  evolving  particle  spectrum  is  determined 
by  solving  the  Fokker-Planck  equation  in  energy  space  for  a 
spatially  homogeneous  ICM,  given  by 


dN(K,  t ) 
dt 


^[W0W0] 

-  ^{  [*„*(*,  t)  -  Atot(K,  t)]N(K,  0} 


N(K,  t) 


Q(K,  t) . 


(21) 


The  quantity  Ajot  represents  the  total  kinetic-energy  loss 
rate,  and  Atot  represents  the  total  energy  gain  rate.  Electrons 
experience  energy  losses  from  synchrotron,  bremsstrahlung, 
Compton,  and  Coulomb  processes.  The  total  continuous 
energy-loss  rate  for  protons  is  due  to  the  effects  of  Coulomb 
processes.  In  addition,  protons  experience  catastrophic 
energy  losses  due  to  proton-proton  collisions  ( i  =  pp), 
proton-7  collisions  ( i  =  p 7),  and  diffusive  escape  (i  =  d)  on 
the  timescale  r,(Ap,f).  The  spectra  of  the  secondary  elec¬ 
trons  are  calculated  from  the  pion-decay  products  (see 
Appendix  A.3),  and  are  subject  to  the  same  energy  losses  as 
the  primary  electrons.  The  calculation  of  the  energy-loss 
rates  accounted  for  by  the  quantity  ATtot  for  both  electrons 
and  protons  is  presented  in  detail  in  Appendix  A.  The  non- 
thermal  energy  spectra  resulting  from  the  temporal  evolu¬ 
tion  of  equation  (21)  will  produce  a  nonthermal  photon 
spectra.  The  method  used  to  calculate  the  synchrotron, 
bremsstrahlung,  Compton,  and  7 r0  7-ray  radiation  compo¬ 
nents  to  the  photon  energy  spectra  is  presented  in  detail  in 
Appendix  B. 

The  integration  of  the  partial  differential  equation  (21) 
poses  a  significant  computational  problem  when  consider¬ 
ing  the  ranges  of  energies  modeled  in  this  paper.  The  prob¬ 
lem  stems  from  the  diffusion  that  occurs  at  low  energies  and 
the  large  energy  losses  that  occur  at  high  energies.  We  found 
that  equation  (21),  when  solved  with  a  Crank-Nicholson 
implicit  finite  differencing  method,  is  unstable  to  an  oscilla¬ 
tion  instability  at  the  high  and  low  energies  because  of  the 
large  values  that  the  convection  derivative  takes  at  both 


1  A  quantitative  prescription  of  the  shock  efficiency  as  a  function  of  the 
injection  index  is  given  by  Miniati  et  al.  (2001);  see  also  Keshet  et  al.  (2003) 
on  this  issue. 
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energy  extremes.  To  solve  equation  (21),  we  settled  on  an 
implicit,  upwind  finite  differencing  method.  This  has  the 
advantage  of  solving  the  oscillation  instability  problem, 
although  with  the  disadvantage  of  a  decrease  in  numerical 
accuracy  for  a  given  step  size. 


2.2. 1 .  Particle  Source  Function 

Nonthermal  particle  distributions  are  injected  into  the 
ICM  according  to  the  function  given  by  equation  (19).  The 
normalization  constant  Q%  is  calculated  as  described  in 
§  2.2.  The  particle  distribution  functions  have  spectral  index 
s(t)  and  exponential  cutoffs  at  kinetic  energy  Ajnax ,  which 
are  determined  by  shock  acceleration  theory  as  described 
below. 

Nonthermal  particles  gain  energy  through  first-order 
Fermi  shock  acceleration  at  the  rate 


Ksh(t)  =  100^^-  MeV  s-1 
P'J 


(22) 


(Stumer  et  al.  1997;  Reynolds  1996),  where  v%(t)  =  v\(t)/ 
(108  cm  s~')  is  the  shock  velocity  in  units  of  103  km  s~',/is 
a  factor  giving  the  mean  free  path  of  a  particle  in  units  of  its 
gyroradius,  and  the  quantity  rj  accounts  for  the  orientation 
of  the  shock  front  relative  to  the  magnetic  field.  In  our 
calculations,  we  set/  =  10  and  /-./  =  1 . 

The  first  constraint  that  we  impose  on  maximum  particle 
energy  is  the  available  time  to  accelerate  particles  to  a 
given  energy.  This  constraint  is  determined  through  the 
expression 

Ain  ax.  t  =  [  dtksh(t)  ,  (23) 

Jo 


where  v\(t)  is  found  from  equation  (4)  and  B  is  assumed  to 
be  constant  throughout  the  cluster  volume.  This  constraint 
is  the  most  important  limitation  on  proton  energy  for  most 
of  the  lifetime  of  the  shock  front. 

Competition  with  energy  losses  sets  a  second  constraint 
on  maximum  particle  energy.  By  equating  the  electron 
energy  gain  rate,  equation  (22),  with  the  sum  of  the  synchro¬ 
tron  and  Compton  energy  losses  that  are  the  dominant 
electron  energy-loss  processes  at  high  energies,  we  obtain  a 
second  energy  constraint  for  electrons,  given  by 


k 


max, 2 


=  2.8  x  10 


,7  "8(0 
Vfrj 


B 


(1  +  z)4  +  10-3Ri- 


1/2 


MeV 


(24) 


where  B  .7  is  the  magnetic  field  in  units  of  10  7  G.  Consider¬ 
ing  synchrotron  losses  only,  the  corresponding  maximum 
energy  for  protons  is 


Aimax,2(syn)  =  2  X  1015 


MO 

\J  B-ifrj 


MeV  . 


(25) 


For  the  parameters  used  in  our  calculations,  proton  energy 
losses  due  to  Compton  and  synchrotron  processes  are 
smaller  than  photomeson  and  photopair  losses  with  the 
CMB.  These  latter  processes  could  in  principle  constrain  the 
maximum  proton  energy  (Berezinskii  &  Grigor’eva  1988) 
but  are  found  to  be  negligible  for  the  parameters  used  in  the 
calculations.  This  is  because  the  timescale  to  accelerate 
protons  to  ultrahigh  energies  (>1019  eV)  requires 
«3  x  10n  [frj/B-jvKt^Kp/lO20  eV)  yr,  which  is  longer 


than  the  Hubble  time  tH  and  much  longer  than  the  timescale 
over  which  the  merger  shocks  persist  in  the  merging  clusters 
(wO.l///).  Although  photomeson  production  does  not  pose 
a  significant  limitation  to  particle  acceleration,  it  is  still 
treated  as  a  catastrophic  loss  in  the  calculation. 

In  order  for  particle  acceleration  to  occur,  the  Larmor 
radius  of  a  particle  must  be  smaller  than  the  size  scale  of  the 
system;  otherwise,  the  particle  escapes  the  system  before 
any  significant  energy  gain  can  occur.  The  third  constraint 
imposed  on  maximum  particle  energy  is  therefore  obtained 
by  requiring  that  the  particle  Larmor  radius  be  less  than  the 
size  scale  of  the  shock  front.  This  constraint  becomes,  for 
both  protons  and  electrons, 


A/ax. 3  —  ' 


eB\  10 
_  _ 

f  ~  f 


13 


-B 


MeV  , 


(26) 


where  A  =  27L  and  e  is  the  electron  charge.  In  all  of  our  cal¬ 
culations,  the  limitations  imposed  by  AjIlaX;i  and  A'max,2  are 
more  restrictive  for  the  electrons.  However,  for  protons, 
Aimax, 3  is  the  strongest  limitation  for  the  highest  energy  pro¬ 
tons  at  the  later  stages  of  the  shock  evolution.  The  maxi¬ 
mum  energy  Kmax(t)  in  equation  (19)  is  the  minimum  of  the 
maximum  energies  calculated  in  equations  (23)-(26).  Note 
that  the  size-scale  limitation  on  maximum  particle  energy 
assumes  that  the  coherence  length  of  the  magnetic  field  is 
larger  than  the  size  scale  of  the  shock. 

The  particle  index  s(t)  is  determined  from  the  collision 
and  sound  speeds  of  the  system.  The  collision  speed  vo(t)  is 
obtained  by  solving  equation  (4)  with  the  initial  radius  of 
the  interaction  given  by  the  sum  of  the  assumed  radii  of  the 
clusters  and  the  initial  collision  velocity  given  by  equation 
(7).  We  assume  that  the  subcluster  is  a  point  mass  and  the 
dominant  cluster  is  described  by  the  isothermal  beta  model 
density  profile  given  in  equation  (3). 

The  shock  compression  ratio  is  calculated  from  the  equa¬ 
tion 


r+  i 

r-i  +  2„/r2(0  ’ 


(27) 


where  M{t)  =  v(t) /cs  is  the  Mach  number  of  the  cluster 
merger  shock  speed,  as  computed  from  equations  (4),  (30), 
and  (31).  The  ratio  of  specific  heats  is  represented  by 
T  =  5/3,  appropriate  to  a  nonrelativistic  monoat omic  ideal 
gas.  In  Figure  2,  the  time  evolution  of  the  compression  ratio 
of  a  shock  front  in  a  merging  cluster  with  masses  M  =  1015 
Mq  and  m  =  1014  M0  is  graphed.  The  radii  of  the  two  sys¬ 
tems  are  assumed  to  be  1.5  and  0.75  Mpc,  and  the  onset  of 
the  merger  is  at  a  redshift  of  z,-  =  0.1.  From  the  shock  com¬ 
pression  ratio,  we  calculate  the  injected  particle  index  5  from 


s(t) 


V{t)  +  2 
V(t)  -  1  ' 


(28) 


The  calculated  shock  speeds  of  cluster  merger  shocks  range 
from  ~1.5  to  3  times  the  ICM  sound  speed,  with  power-law 
indexes  ranging  from  ~5  to  ~2.3.  Figure  2  also  shows  the 
evolution  of  the  typical  power-law  index  as  a  function  of 
time. 

In  the  merger  event  between  two  clusters  of  galaxies,  as 
described  in  §  2.1,  two  shock  fronts  develop.  The  shock 
front  propagating  into  the  dominant  cluster,  which  we  call 
the  forward  shock,  has  Mach  number  J4\[t).  Similarly,  the 
shock  front  propagating  into  the  merging  cluster  has  Mach 
number  Jljif)-  Writing  down  the  compression  ratio  for  the 
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Fig.  2. — Compression  ratio  of  the  forward  ( dashed  curve)  and  reverse 
(solid  curve)  shocks  as  calculated  by  the  method  described  in  §  2.2.1.  The 
particle  injection  spectral  index  for  the  forward  (dotted  curve )  and  reverse 
(dot-dashed  curve)  shocks  is  associated  with  the  compression  ratio  and  is 
calculated  for  a  merger  system  consisting  of  a  dominant  cluster  of  mass 
1015  Me  and  a  poor  subcluster  of  mass  1014  M0.  Initial  onset  of  the  shock 
begins  at  a  radius  of  1 .5  Mpc  from  the  cluster  center  with  an  initial  redshift 
of  z,-  =  0.1.  Particle  acceleration  occurs  over  a  period  facc  =1.1  Gyr  for 
these  parameters,  after  which  injection  stops.  Acceleration  at  the  reverse 
shock  ends  after  the  reverse  shock  passes  through  the  merging  cluster. 


forward  \%>\  (?)]  and  reverse  [W 2{t )]  shocks  from  equation 
(27),  we  get 


<?i(0  = 
V2(t)  = 


»i(Q  =  4 

v  ~  3[1  -J/\2(t)\  ’ 
i/2(0  4 

~T  =  3[l-.#2-2(0]  ’ 


(29) 


where  u\{t)  and  u2{t)  are  the  velocities  of  the  forward  and 
reverse  shocks,  respectively.  In  addition,  v  and  v  are  the 
velocities  of  the  shocked  fluid  for  the  forward  and  reverse 
shocks,  respectively.  Noting  that  v  =  vo(t)  -  v,  where  Vo  is 
the  infall  velocity  of  the  merging  cluster,  and  solving  for  the 
Mach  number,  we  get 


2  v0(t)  -  v 

3  c2 


[pO(0  -  v] 


(30) 


where  c\  is  the  sound  speed  in  the  dominant  cluster,  c2  is  the 
sound  speed  in  the  merging  cluster,  and  the  Mach  number  is 
defined  to  be  ,2/t'i  ,2-  The  values  of  <7  and  C2  are 

calculated  from  equation  (8)  with  temperatures  given  by 
equation  (9)  for  clusters  of  mass  M\  and  M2.  The  value  of  v 
is  calculated  by  solving 

/q«i(f)  1  +  3 

R2n2(t)  1  +  3  Jt22{i) 


Vo (t)  -  V 


(31) 


where  n\(t)  and  n2{t)  are  the  ICM  number  densities  in  the 
dominant  cluster  and  merging  cluster  at  the  positions  of  the 
forward  and  reverse  shocks,  respectively,  and  are  calculated 
according  to  equation  (3).  We  take  the  mean  mass  per  par¬ 
ticle  /q  =  /z2  =  0.6 mp  in  the  two  clusters.  Equation  (31)  was 
derived  from  the  shock  jump  conditions  by  equating  the 
energy  densities  of  the  shocked  fluids  at  the  contact  dis¬ 
continuity.  The  compression  ratios  are  then  calculated  from 
equation  (27).  These  values  are  plotted  in  Figure  2. 


3.  RESULTS 

We  present  results  of  the  numerical  model  described  in  §  2 
to  simulate  the  temporal  evolution  of  nonthermal  particles 
accelerated  by  shocks  formed  in  merging  clusters  of  gal¬ 
axies.  From  these  particle  distributions,  we  are  able  to  inves¬ 
tigate  the  expected  nonthermal  photon  spectra  for  energies 
ranging  from  ^lO-9  to  ~1019  eV  that  cover  astrophysically 
interesting  frequencies.  At  photon  energies  <4  x  1 0  3  eV 
photons  («10  MHz),  ground-based  radio  observations  are 
subject  to  the  ionospheric  absorption  cutoff.  At  photon 
energies  >  1012  eV  (depending  precisely  on  source  redshift), 
the  universe  becomes  opaque  to  7-rays  due  to  pair  produc¬ 
tion  attenuation  on  the  diffuse  infrared  radiation  fields.  The 
electron-positron  pairs  produced  by  the  ^1014  to  ~1019  eV 
photons  provide  a  secondary  source  of  synchrotron  and 
Compton-scattered  radiation  on  the  Mpc  galaxy-cluster 
scale,  but  this  cascade  radiation  is  found  to  provide  only  a 
minor  nonthermal  synchrotron  component  to  the  nonther¬ 
mal  UV/soft  X-ray  emission  and  is  not  included.  To 
calculate  detectability  of  clusters  with  ground-based 
Cerenkov  telescopes,  77  absorption  of  the  high-energy  7- 
rays  must,  however,  be  considered. 

The  energy-loss  timescales  of  nonthermal  electrons  for 
various  parameters  assumed  in  cluster  environments  are 
shown  in  Figure  3.  At  low  energies,  Coulomb  losses  domi¬ 
nate.  A  maximum  in  the  energy-loss  timescale  is  found  in 
the  energy  range  from  « I  to  several  hundred  MeV.  At 
higher  energies,  Compton  and  synchrotron  losses  dominate. 
At  the  highest  energies,  the  Compton  loss  rate  declines 
because  of  the  onset  of  the  Klein-Nishina  reduction  in  the 
Compton  cross  section,  and  synchrotron  losses  become  the 
dominant  loss  mechanism  (see  also  Sarazin  1999).  The  elec¬ 
tron  loss-rate  from  double  Compton  scattering  is  always 
much  smaller  than  synchrotron  losses  in  our  study  and  can 
be  safely  neglected  (see  §  A.  1 ). 

Figure  4  shows  the  nonthermal  particle  energy  spectra  at 
different  times  during  the  merger  of  a  dominant  cluster  of 
mass  M\  =  1015  M0  and  radius  Ri  =  1.5  Mpc  with  a  sub¬ 
cluster  of  mass  M2  =  1014  M0  and  a  radius  R2  =  0.75  Mpc, 
respectively.  The  quantities  M \ ,  M2,  Ru  and  R2  are  defined 
in  §  2. 1 . 1 .  The  assumed  core  radius  rc  for  the  dominant  clus¬ 
ter  is  250  kpc.  The  ICM  density  of  the  dominant  cluster  is 
assumed  to  follow  equation  (3)  with  a  central  number  den¬ 
sity  no  =  10-3  cm-3  and  (3  =  0.75.  These  parameters  are 
consistent  with  the  Lx-7x  relation  (eq.  [12]).  A  uniform 
magnetic  field  of  strength  1  /uG  was  assumed  throughout 
the  cluster.  The  galaxy  cluster  environment  and  its  evolu¬ 
tion  is  calculated  using  the  cosmological  parameters 
(fl0,  flRl  flA)  =  (0.3,  0.0,  0.7).  All  models  assume  a  Hubble 
constant  Ho  =  70  km  s-1  Mpc-1.  We  also  assume  a 
hydrogen-to-helium  number  ratio  of  10. 

Irrespective  of  the  specific  parameters  of  the  simulation, 
we  see  that  at  energies  <  1  MeV,  the  electron  and  proton 
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Fig.  3. — Instantaneous  energy  loss  timescale  for  electrons  in  a  merging  cluster  scenario  at  redshift  z  =  0.0.  In  the  left  panel,  the  proton  density 
n  =  10' 3  cm'3,  and  B  =  1  // G.  Only  the  dominant  CMB  radiation  field  is  considered.  In  the  right  panel,  n  =  10'6  cm'3,  and  5  =  0.1  /rG. 


spectra  display  low-energy  cutoff's  due  to  Coulomb  diffusion 
and  energy  losses  (the  thermal  particle  distributions  are  not 
shown).  Electrons  with  energies  greater  than  ~1  GeV  lose  a 
significant  fraction  of  their  energy  on  Gyr  timescales.  The 
dominant  energy-loss  mechanism  for  these  electrons  is 
Compton  scattering,  as  Thomson  losses  on  the  CMB  domi¬ 
nate  synchrotron  losses  when  B-4  3(1  +  z)2  /iG  (see  Fig.  3). 
The  dominant  energy-loss  mechanism  for  protons  is  secon¬ 
dary  nuclear  production,  which  operates  on  a  timescale  of 
about  2 tH  for  an  ambient  density  of  10~3  cm-3.  Conse¬ 
quently,  the  high-energy  cutoff  in  the  proton  distribution 
function  results  from  the  available  time  to  accelerate  par¬ 
ticles  and,  at  the  longest  times,  competition  of  particle  accel¬ 
eration  with  catastrophic  pp  losses.  Even  though  particles 
are  injected  with  a  range  of  power-law  indexes,  most  of  the 
injection  power  occurs  when  the  spectral  index  is  hardest, 
yielding  a  proton  spectrum  that  can  be  accurately  described 
by  a  single  power  law  (compare  Fig.  2)  with  a  number  index 
of  «— 2.5  for  the  forward  shock  and  «— 2.2  for  the  reverse 
shock.  The  electron  bremsstrahlung  and  secondary  nuclear 
production  assume  a  target  proton  density  of  1(F3  cm~3. 


Separate  radiation  components  produced  by  nonthermal 
particles  are  shown  in  Figure  5  for  a  shock  formed  in  the 
merger  of  1014  and  1015  M0  clusters,  with  the  interaction 
beginning  at  zt  =  0.3.  The  system  has  been  evolved  for  0.8 
Gyr  (z  =  0.22)  in  Figures  5 a  and  5b  and  for  3.42  Gyr  (z  =  0) 
in  Figure  5c.  Figures  5 a  and  5b  correspond  to  times  before 
particle  injection  has  stopped.  The  cluster  magnetic  field 
strength  B  =  1.0  pG  in  Figures  5a  and  5c  and  B  =  0.1  pG  in 
Figure  5b.  The  radiation  components  of  the  primary  and 
secondary  electrons  are  shown  by  the  short-dashed  and  dot¬ 
ted  curves,  respectively,  and  the  7r°-decay  7-rays  are  shown 
by  the  dotted  curves.  The  curves  are  labeled  by  the  various 
radiation  processes,  and  the  total  spectral  energy  distribu¬ 
tions  are  shown  by  the  solid  curves.  The  synchrotron  and 
Compton  radiations  from  primary  electrons  dominate  the 
emission  at  times  before  primary  injection  stops,  although 
7r°-decay  7-rays  make  a  dominant  contribution  above  100 
TeV.  At  late  times  after  particle  injection  ceases,  the  cooling 
primary  electrons  provide  the  dominant  contribution  at  the 
lowest  radio  frequencies  and  at  optical/UV  frequencies,  but 
secondary  electrons  can  make  the  dominant  contribution  at 


Fig.  4. — Primary  electron  (left)  and  proton  (right)  energy  spectra  for  a  shock  formed  by  a  merging  cluster  using  a  standard  parameter  set  with  M\  =  1015 
M0,  M2  =  1014  M0,  and  B  =  1.0  fj.G.  The  interaction  begins  at  redshift  z,-  =  0.3  and  is  evolved  to  the  present  time  (z  =  0.0;  solid  curve )  with  spectra  shown  at 
the  various  times  given  in  the  figure.  The  redshifts  corresponding  to  the  time  delays  of  0.1,  0.5,  0.8,  1.0,  1.2,  1.5,  and  3.42  Gyr  arez  =  0.289, 0.248, 0.218, 0.199, 
0.180,  0.153,  and  0.0,  respectively.  Particle  injection  ends  at  /acc  =  1.1  Gyr  (see  discussion  following  eq.  [19]). 
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Fig.  5. — Separate  radiation  components  produced  by  nonthermal  particles  accelerated  by  a  shock  initiated  at  z,-  =  0.3  resulting  from  the  merger  between 
1014  and  1015  M0  clusters.  The  separate  radiation  components  are  labeled  in  the  figures.  Emission  from  the  primary  electrons  are  represented  by  the  dashed 
curves,  and  emission  from  secondary  electrons  and  positrons  and  7r°  7-rays  are  represented  by  the  dotted  curves.  The  solid  curves  are  the  total  nonthermal 
photon  spectra.  Panels  a  and  b  show  the  nonthermal  spectrum  at  0.8  Gyr  (z  =  0.22),  and  panel  c  at  3.42  Gyr  (z  =  0.0)  after  the  merger  shock  forms.  Magnetic 
field  strengths  are  B  =  1.0  /iG  for  panels  a  and  c  and  B  =  0.1  fiG  for  panel  b.  Panels  a  and  b  correspond  to  times  before  the  termination  of  the  nonthermal 
particle  injection.  Panel  c  is  approximately  2.3  Gyr  after  particle  injection  has  ceased. 


1.4  GHz.  At  7-ray  energies,  primary  bremsstrahlung  and 
secondary  7r°-decay  7-rays  make  the  dominant  contribution 
at  late  times.  Note  that  structure  in  some  of  the  emission 
components  is  due  either  to  knock-on  electrons  or  to  the 
separate  effects  of  the  stellar,  X-ray,  and  CMB  Compton 
energy  losses  on  the  electron  energy  spectra.  As  can  be  seen, 
the  intensity  of  the  synchrotron  radio  emission  is  very 
sensitive  to  the  magnetic  field  strength. 

Figure  6  shows  the  temporal  evolution  of  the  total  photon 
spectra  radiated  by  a  shock  formed  in  the  merger  of  1014  and 
1015  Mq  clusters  with  B  =  1.0  juG.  The  photon  spectra  are 
shown  at  various  times  for  shocks  that  begin  at  redshifts 
Zj  =  0.03,  0.05,  0.1,  and  0.3.  Each  shock  is  evolved  to  the 
present  time  (z  =  0.0),  and  the  radiation  spectrum  at  the 
present  time  is  represented  by  the  solid  curve.  Typical  life¬ 
times  for  the  shocks  are  ~1  Gyr.  As  a  result,  only  shock 
fronts  with  z,  >  0.1  show  late-time  evolution  of  the  total 
photon  spectra. 

At  early  times  well  before  particle  injection  ceases  (<0.5 
Gyr),  the  luminosity  increases  as  the  shock  front  accelerates 
particles  and  the  particles  accumulate.  During  this  epoch, 
the  injection  spectrum  is  very  soft,  and  nonthermal  brems¬ 


strahlung  and  7r°-radiation  can  dominate  the  production  of 
photons  at  high  energies  because  of  the  relatively  few  high- 
energy  electrons.  The  peaks  in  the  spectra  at  t  <  0.3  Gyr 
spectra  are,  from  low  to  high  energies,  due  to  Compton, 
bremsstrahlung,  and  7r°-radiation.  The  enhancement  at  low 
frequencies  is  due  to  nonthermal  synchrotron  radiation.  As 
the  front  ages  and  the  electron  spectrum  hardens  (0.5-1 
Gyr),  the  number  of  high-energy  electrons  increases,  and 
the  Compton  process  strongly  dominates  the  nonthermal 
emission  at  infrared  energies  and  above. 

At  times  between  tcou  and  tacc  (between  ~0.9  and  1.1  Gyr 
in  Fig.  6),  the  shock  front  propagates  to  the  outer  regions  of 
the  clusters,  and  the  particle  injection  rate  decreases.  The 
number  of  high-energy  electrons  begins  to  decrease  as  their 
previous  numbers  are  no  longer  replenished.  The  conse¬ 
quence  of  this  behavior  is  a  softening  in  the  synchrotron 
and  high-energy  Compton  emission.  This  is  seen  clearly  in 
the  Zi  =  0.1  and  r,  =  0.3  panels.  As  the  Compton  emission 
declines,  the  7-rays  from  the  p  +  p  — >  n0  — >  27  process 
become  increasingly  evident  at  times  >1.1  Gyr. 

The  shock  front  terminates  at  ta cc  after  propagating 
through  the  cluster.  At  these  late  times,  particle  injection 
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Fig.  6. — Temporal  evolution  of  the  total  spectral  energy  distribution  for  a  cluster  merger  shock  that  begins  at  redshift  z,  =  0.03  (top  left),  z,-  =  0.05  ( top 
right),  z,-  =  0.1  (bottom  left),  or  z,  =  0.3  (bottom  right).  Each  shock  is  evolved  to  the  present  (solid  curve).  Intermediate  times  are  indicated  in  each  panel.  In  all 
cases  the  shock  survives  for  ~  1 . 1  Gyr.  Cluster  merger  shocks  that  develop  at  redshifts  >0.1  do  not  persist  to  the  present  epoch.  The  thin  solid  curve  shows  the 
expected  thermal  bremsstrahlung  for  a  1015  MQ  cluster,  whose  temperature  is  given  by  eq.  (9).  The  figures  with  two  thin  solid  curves  represent  the  thermal 
bremsstrahlung  at  z,  (lower  curve)  andz  =  0.0  (higher  curve).  The  integrated  luminosity  of  the  thermal  bremsstrahlung  is  «1045  ergs  s-1. 


ceases,  and  the  remaining  population  of  nonthermal 
particles  evolve  according  to  the  energy  loss  and  diffusion 
processes  described  in  §§  A1  and  A2,  with  no  freshly  acceler¬ 
ated  electrons  to  repopulate  the  particle  distribution.  The 
system  quickly  becomes  devoid  of  primary  electrons,  so  that 
the  bremsstrahlung  and  pion  emission  dominate  at  HXR 
energies  and  above.  The  most  extreme  case  is  shown  in  the 
zt  =  0.3  panel  of  Figure  6.  The  spectrum  at  the  present 
epoch  shows  three  maxima  produced  by  synchrotron  radia¬ 
tion  at  kHz  frequencies,  Compton  scattering  near  100  eV, 
and  combined  nonthermal  bremsstrahlung  and  7ro-decay 
emission  near  100  MeV.  The  plateau  in  the  radio  emission 
at  late  times  is  formed  by  synchrotron  emission  from 
secondary  electrons  (1.4  GHz  corresponds  to  6  x  10  6  eV). 

Also  shown  in  Figure  6  is  the  thermal  bremsstrahlung 
emission  from  the  ICM  using  the  form  for  the  emission 
spectrum  given  by  Rybicki  &  Lightman  (1979).  The  temper¬ 
ature  and  the  luminosity  of  the  thermal  bremsstrahlung  is 
computed  from  equations  in  §  2.1.2.  Thermal  bremsstrah¬ 
lung  generally  dominates  the  emission  at  10  eV  to  10  keV 
energies.  It  should  be  noted  that  the  relations  for  tempera¬ 
ture  and  luminosity  of  the  thermal  bremsstrahlung  compo¬ 
nent  do  not  take  into  consideration  the  merging  process  of 


the  system.  Computer  simulations  of  expected  X-ray 
temperatures  and  luminosities  of  merging  clusters  indicate 
an  increase  in  the  thermal  bremsstrahlung  luminosity  and 
temperature  (Ricker  1998;  Schindler  &  Muller  1993).  There¬ 
fore,  the  thermal  bremsstrahlung  curves  presented  in  Fig¬ 
ures  6  and  7  very  likely  underestimate  the  true  X-ray 
luminosity  and  temperature  of  the  gas  during  the  merger 
event.  At  times  ;>>  tacc,  the  relation  between  the  initial  tem¬ 
perature  and  luminosity  will  depend  on  mass  ratio  of  the 
merger  components,  the  dark-matter  density  profile  of  the 
clusters,  and  impact  parameter  of  the  merger  event  (Ricker 
&  Sarazin  2001 ). 

Figure  7  shows  a  calculation  for  the  same  system  as  in 
Figure  6,  except  for  a  cluster  with  a  magnetic  field  of  0. 1  /uG. 
We  consider  these  two  values  because  of  the  uncertainty 
regarding  the  mean  magnetic  fields  in  clusters  of  galaxies. 
Despite  Faraday  rotation  measurements  indicating  mag¬ 
netic  fields  strengths  B>  1 .0  juG  (Clarke,  Kronberg,  & 
Bohringer  2001),  a  comparison  of  radio  and  Compton 
luminosities  suggests  field  strengths  B«0.1  /iG  (Fusco- 
Femiano  et  al.  1999).  Possible  explanations  for  this 
discrepancy  may  arise  from  an  observational  biasing  of  the 
Faraday  rotation  measurement  of  the  magnetic  fields — 
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Fig.  7. — Same  as  Fig.  6,  except  with  an  assumed  cluster  magnetic  field  of  0. 1  /iG 


weaker  magnetic  fields  mean  brighter  Compton  emission 
for  a  fixed  radio  luminosity — and  the  possibility  that  the  dis¬ 
tributions  of  the  magnetic  field  and  nonthermal  particles  are 
anticorrelated  due  to  the  accumulation  of  electrons  in  weak 
magnetic  field  regions  (Petrosian  2001).  The  major  change 
between  the  two  figures  is  an  increase  in  the  synchrotron 
radio  emission,  which  is  roughly  proportional  to  the  square 
of  the  magnetic  field.  The  system  with  the  lower  magnetic 
field  has  a  slightly  brighter  Compton  component. 

Light  curves  at  various  observing  frequencies  are  shown 
in  Figure  8  using  the  parameters  for  our  standard  case  given 
by  Figure  6  with  5=1.0  //G  and  a  cluster  merger  beginning 
at  Zi  =  0.3.  The  light  curves  of  the  nonthermal  radiation 
exhibit  a  common  behavior  independent  of  frequency.  At 
early  times,  the  spectral  power  rises  rapidly  as  the  clusters 
merge.  The  peak  emission  occurs  when  the  centers  of  mass 
of  the  two  clusters  pass  at  Coin  after  which  the  emission 
exhibits  a  slow  decay  and  approaches  a  plateau  at  times 
f>4cc  when  particle  injection  has  stopped.  The  rate  of 
decay  of  the  emission  increases  with  radio  frequency  due  to 
the  stronger  cooling  of  the  higher  energy  electrons,  so  that 
the  decay  is  slowest  at  lower  frequencies.  Synchrotron  emis¬ 
sion  from  secondary  electrons  forms  the  late-time  plateaus 
at  radio  energies.  This  behavior  is  also  apparent  for  the 
HXR  emission,  although  it  is  formed  by  primary  brems- 
strahlung  and  both  primary  and  secondary  Compton 


radiation  at  late  times.  At  7-ray  energies,  the  7r°-decay  emis¬ 
sion  forms  a  plateau  of  emission  that  dominates  soon  after 
t  =  Lee  ■ 

In  Figure  9  we  show  the  spectral  energy  distributions  and 
spectral  components  for  a  case  where  the  efficiency 
ijp  =  10%  for  protons  and  rje  =  1%  for  electrons.  An 
enhancement  of  the  number  flux  of  nonthermal  protons 
compared  to  electrons  is  expected  in  nonrelativistic  shock 
acceleration  because  protons,  unlike  electrons,  resonate 
with  Alfven  waves  at  all  energies  (Baring  et  al.  1999),  and 
the  Larmor  radii  of  protons  in  the  tail  of  the  Maxwellian 
distribution  of  the  background  thermal  plasma  is  greater 
than  that  of  the  electrons,  so  that  protons  can  more  effec¬ 
tively  scatter  across  the  shock  front.  This  theoretical  expect¬ 
ation  is  borne  out  by  the  large  proton-to-electron  ratio  in 
the  cosmic  rays.  The  enhanced  proton  contribution  produ¬ 
ces  structure  in  the  7-ray  spectrum  even  at  t  w  ?coU,  when  the 
emission  is  brightest.  The  7r°-decay  feature  will  be  clearly 
evident  in  the  spectrum  of  merging  clusters  of  galaxies  for 
such  a  large  proton-to-electron  efficiency  ratio. 

Calculations  of  the  hardest  particle  injection  spectral 
index  .vmm  formed  in  cluster  merger  shocks  are  shown  in 
Figure  10  as  a  function  of  the  larger  mass  M\  of  the  two 
clusters,  with  the  subcluster  mass  M2  =  1014  M0.  To  derive 
these  results,  we  used  equations  in  §  2.1.1  to  calculate  the 
dynamics  of  the  clusters,  equations  in  §  2. 1 .2  to  calculate  the 
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Fig.  8. — Light  curves  at  various  observing  frequencies  produced  by  a  shock  formed  in  a  merger  between  1014  and  1015  M0  clusters  that  begins  at  r,  =  0.3 
and  is  evolved  to  the  present  epoch  (t  =  3.42  Gyr).  All  light  curves  are  for  a  magnetic  field  strength  of  B  =  1.0  /rG  unless  otherwise  noted.  Radio  light  curves  in 
janskys  are  given  at  15  MHz,  30  MHz,  74  MHz,  1 10  MHz,  and  1.4  GHz  in  the  left  panel,  and  light  curves  in  energy  flux  units  are  given  at  40-80  keV, 
>100  MeV,  >1  GeV,  and  >100GeV  in  the  right  panel.  The  15  MHz  light  curve  is  also  calculated  with  a  magnetic  field  strength  of  B  =  0.1  //G  (dotted  curve). 


sound  speed,  and  equations  (27)  ( 3 1 )  to  calculate  the  spec¬ 
tral  index  from  the  Mach  number.  We  also  assume  that  the 
onset  of  the  merger  begins  at  redshift  z-L  =  0.1;  softer  injec¬ 
tion  indexes  are  obtained  for  mergers  at  larger  values  of  z,- 
because  of  the  smaller  maximum  separations  d  at  earlier 
times,  given  by  equation  (1).  We  calculate  ,smm  for  various 
values  of  rc  and  /3.  The  values  (rc,/3)  =  (0.05,  0.8),  (0.05, 
0.45),  (0.5,  0.8),  and  (0.5,  0.45)  roughly  correspond  to  the 
extrema  in  the  range  of  these  parameters  measured  for  45 
X-ray  clusters  observed  with  ROSAT  (Wu  &  Xue  2000). 
Also  shown  are  values  of  ,smjn  for  (rc,/3)  =  (0.179,  0.619), 
which  are  the  average  values  of  these  parameters  for  the  45 
X-ray  clusters,  and  (rc,f 3)  =  (0.25,  0.75),  which  are  the 
standard  parameters  used  in  the  calculations. 

The  uppermost  curves  with  large  />- values  display  a 
maximum.  At  larger  values  of  M\ ,  .vmm  hardens  because  of 
the  stronger  gravitational  potential.  Although  5'mm  spans  a 
wide  range  of  values,  the  value  of  ,vmm  corresponding  to  the 


average  values  of  rc  and  b  lies  in  the  range  2.1  <Amin<2.4, 
and  more  realistically  in  the  range  2.2<Amjn <2.4  because 
clusters  with  masses  3  x  1015  to  1016  M0  are  very  rare. 
Because  z,-  =  0.1,  there  is  a  constraint  on  the  available  time 
for  the  interaction  to  occur.  At  low  values  of  M\ ,  the  gravi¬ 
tational  interaction  is  small  so  that  the  interaction  takes  a 
longer  time  to  occur.  The  features  in  the  curves  between 
1014  and  1014  5  Me  are  a  result  of  this  effect. 

4.  DISCUSSION 

We  have  presented  results  of  a  computer  simulation 
designed  to  calculate  the  temporal  evolution  of  nonthermal 
particles  accelerated  by  a  shock  front  formed  between  two 
merging  clusters  of  galaxies.  We  then  calculated  the 
expected  total  nonthermal  photon  spectra  resulting  from 
synchrotron,  bremsstrahlung,  and  Compton  radiation 
from  both  primary  and  secondary  electrons  and  7ro-decay 


Fig.  9. — Same  as  the  z,-  =  0.3  case  in  Fig.  6,  except  that  r/,,  =  1%  and  r)p  =  10%.  Right  panel  shows  the  separate  radiation  components  at  time  t  =  0.8  Gyr 
using  the  same  notation  as  Fig.  5,  and  the  left  panel  shows  the  total  spectral  power  at  various  times. 
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Fig.  10. — Calculations  of  the  hardest  particle  injection  spectral  indexes 
■Smin  formed  in  cluster  merger  shocks  as  a  function  of  the  larger  mass  M\  of 
the  two  clusters,  for  various  values  of  rc  and  /3.  The  values  of  rc  are  given  in 
Mpc.  The  minimum  spectral  index  for  the  forward  and  reverse  shock  is 
shown  in  the  figure.  The  mass  of  the  subcluster  is  Mi  =  1O14M0. 


7-radiation  from  accelerated  protons.  Our  results  apply  to 
the  shocks  formed  in  the  merger  of  two  clusters  of  galaxies 
and  not  to  the  accretion  shocks  formed  by  spherical  infall 
of  matter  during  the  process  of  structure  formation. 
Theoretical  studies  suggest  that  accretion  shocks  should 
exist  (Bertschinger  1985;  Miniati  et  al.  2000)  and  that  each 
cluster  itself  will  be  the  center  of  its  own  accretion  flow. 
Although  X-ray  and  radio  observations  (Kassim  et  al.  2001; 
Markevitch  &  Vikhlinin  2001;  Fusco-Femiano  et  al.  2003b; 
Bagchi  2003;  Valtchanov  et  al.  2002)  of  halos  and  relics  pro¬ 
vide  strong  evidence  for  the  existence  of  merger  shocks, 
there  is  also  observational  evidence  for  accretion  shocks 
(Ensslin  et  al.  1998, 2001;  Slee  et  al.  2001;  Ensslin  &  Briiggen 
2002).  We  have  focused  on  merger  shocks  in  this  study.  Pre¬ 
vious  work  by  Fujita  &  Sarazin  (2001)  would  classify  our 
simulations  as  accretion  events,  and  they  showed  the  lumi¬ 
nosity  of  these  accretion  shocks  to  be  ~  1 043  ergs  s  1  when 
Mi  =  1  x  1015  M0,  which  is  consistent  with  the  results  of 
our  simulations. 

Beppo-SAX  and  Rossi  X-Ray  Timing  Explorer  observa¬ 
tions  show  the  presence  of  a  HXR  excess  at  energies  >  20 
keV  from  A1656  (Coma)  (Fusco-Femiano  et  al.  1999; 
Rephaeli  et  al.  1999),  A2256  (Fusco-Femiano  et  al.  2000), 
and  A3667  (Fusco-Femiano  et  al.  2003a).  Our  simulations 
support  a  Compton-scattering  CMB  interpretation  of  the 
HXR  (Rephaeli  1979)  and  EUV  (Hwang  1997;  Ensslin  & 
Biermann  1998;  Lieu  et  al.  1999b;  Atoyan  &  Volk  2000) 
emission.  As  shown  in  Figures  6  and  7,  the  nonthermal  radi¬ 
ation  from  the  merger  shock  during  the  postcollision  evolu¬ 
tion  (>tCoii)  shows  an  excess  above  the  thermal 
bremsstrahlung  in  both  the  HXR  and  EUV  regions  of  the 
electromagnetic  spectrum.  For  the  standard  parameters 
used,  the  calculated  luminosity  of  the  nonthermal  compo¬ 
nent  is  Lnt  ~  1043  for  the  HXR  emission  between  the  20 


and  80  keV  energies.  For  the  Coma  Cluster,  BeppoSAX 
observations  indicate  that  the  HXR  flux  is  2.2  x 
10~n  ergs  cm-2  s_1,  which  corresponds  to  an  observed 
luminosity  of  ~1043ergs  s  1  in  the  20-80  keV.  Our  simula¬ 
tions  for  the  merger  model  between  the  times  /con  and  ?acc 
are  consistent  with  this  observed  value.  Our  model  does 
match  well  the  calculated  mass  ratio  of  the  Coma  Cluster 
merger  model  (Colless  &  Dunn  1996).  In  A2256,  the 
observed  masses  of  the  merging  clusters  are  1.6  x  1015  and 
5.1  x  1014  M0  (Berrington  et  al.  2002).  The  observed  HXR 
flux  is  1.2  x  10~n  ergs  cm-2  s_1,  which  corresponds  to  an 
HXR  luminosity  ~  1 044  ergs  s  1  (Fusco-Femiano  et  al. 
2000).  By  scaling  our  models  to  the  observed  masses  of  the 
merging  clusters,  the  expected  HXR  luminosity  is  ~1044 
ergs  s  '1,  which  is  consistent  with  the  observed  values  seen  in 
A2256. 

The  Extreme  Ultraviolet  Explorer  ( EUVE )  has  detected 
EUV  emission  in  the  ~60  250  eV  band  for  a  number  of  clus¬ 
ters.  These  clusters  include  Virgo,  Coma,  Fornax,  A2199, 
A1795,  and  A4059  (Lieu  et  al.  1996a,  1996b,  1999a,  1999b; 
Berghofer,  Bowyer,  &  Korpela  2000b;  Bowyer,  Lampton,  & 
Lieu  1996;  Bowyer  &  Berghofer  1998;  Mittaz  et  al.  1998; 
Bonamente,  Lieu,  &  Mittaz  2001;  Kaastra  et  al.  1999). 
Note,  however,  that  no  detection  of  EUV  emission,  other 
than  that  which  could  be  attributed  to  thermal  tail  X-ray 
emission,  is  reported  from  analysis  of  EUVE  observations 
for  Fornax  (Bowyer,  Korpela,  &  Berghofer  2001),  A2199 
and  A1795  (Berghofer  &  Bowyer  2002),  and  A4059 
(Berghofer,  Bowyer,  &  Korpela  2000a).  The  measured  EUV 
excess  in  the  60-250  eV  band  for  Coma  and  A 1795  is 
2  x  1043  ergs  s_1  (Lieu  et  al.  1996a)  and  3  x  1043  ergs  s_I, 
respectively  (Bonamente  et  al.  2001).  The  calculated  EUV 
luminosity  of  the  nonthermal  component  from  our  simula¬ 
tions  is  L\t  «  1043  for  these  EUV  energies  and  is  consistent 
with  both  clusters.  As  mentioned  before,  our  simulations 
are  a  good  match  to  the  Coma  Cluster  merger  scenario.  The 
masses  of  the  merging  clusters  in  A1795  are  uncertain 
(Oegerle  &  Hill  1994),  and  confirming  the  results  of  our  sim¬ 
ulations  with  the  observed  luminosity  of  the  EUV  excess  of 
A 1795  is  difficult.  An  EUV  excess  of  luminosity  ~6  x  1042 
ergs  s  1  was  detected  from  the  cluster  A2199  (Kaastra  et  al. 
1999).  However,  A2199  lacks  strong  evidence  for  a  current 
merger  event.  If  any  mergers  occurred  in  the  formation 
history  of  A2199,  they  occurred  long  enough  in  the  past  to 
allow  the  cluster  to  revirialize.  It  is  possible  that  an  adjacent 
cluster  A2197  is  beginning  to  interact  with  A2199  (Oegerle 
&  Hill  2001),  but  seems  unlikely  that  it  is  the  origin  of  the 
EUV  excess  in  the  central  400  kpc  of  the  cluster.  Depending 
on  the  masses  of  the  merging  clusters,  the  decreased  lumi¬ 
nosity  of  the  EUV  emission  in  A2199  is  consistent  with  a 
cooling  nonthermal  particle  population.  Radio  and  X-ray 
observations  of  the  central  regions  of  A2199  indicate  it  is 
complex  (Owen  &  Eilek  1998),  and  it  is  possible  other 
mechanisms  are  producing  the  EUV  excess  emission. 

We  treated  in  detail  the  merger  of  a  1014  M0  cluster  with  a 
1015  M0  cluster  having  values  of  rc  and  (3  in  equation  (3) 
approximately  equal  to  the  mean  values  found  in  a  survey 
of  ROSAT  clusters  of  galaxies  (Wu  &  Xue  2000).  We  found 
that  the  hardest  spectral  index  reached  by  the  merger  shock 
in  this  system  is  .vmm  «  2.2  for  zt  =  0.1.  Because  rc  and  A  are 
not  correlated  with  cluster  temperature  7x  and  therefore 
cluster  mass,  as  can  easily  be  seen  from  the  results  of  Wu  & 
Xue  (2000),  our  model  results  are  reasonably  representative 
of  typical  observed  cluster  mergers. 
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We  consider  the  detectability  of  nonthermal  synchrotron 
emission  with  the  Low  Frequency  Array  (LOFAR).2 
LOFAR  is  a  multielement  interferometric  low-frequency 
radio  telescope  designed  for  high-resolution  radio  imaging 
in  the  10-300  MHz  frequency  range.  Figure  8  includes  sev¬ 
eral  frequencies  observable  by  LOFAR  showing  the  optimi¬ 
zation  of  the  LOFAR  telescope  for  detecting  these  objects. 
Estimated  sensitivities  for  the  LOFAR  telescope  range  from 
1.6  mJy  at  15  MHz  to  40  //Jy  at  1 10  MHz.  The  quoted  sensi¬ 
tivities  are  for  1  hr  exposures  for  a  single  beam  for  a  total  of 
13,365  dual-polarization  dipoles  optimized  in  the  10-90 
MHz  and  213,840  dual-polarization  dipoles  optimized  for 
the  110-220  MHz,  each  with  a  maximum  baseline  of  400 
km.  With  station  diameters  of  65  m,  the  resolution  of 
LOFAR  is  approximately  12"  at  15  MHz  and  2'.'4  at  110 
MHz  at  the  full  extent  of  the  array.  Because  multielement 
interferometric  radio  telescopes  have  trouble  detecting 
large,  extended  sources,  we  must  consider  the  angular  size 
of  the  shock  front.  Given  a  size  scale  of  ^1  Mpc,  these  fea¬ 
tures  will  easily  be  resolved  at  all  frequencies  visible  to 
LOFAR  at  a  distance  of  100  Mpc.  As  a  worst-case  scenario, 
we  consider  the  extent  of  the  array  that  is  only  capable  of 
observing  the  shock  front  as  a  point  object.  For  all  but  the 
most  distant  shock  fronts  and  the  lowest  frequencies,  the 
central  2  km  of  the  LOFAR  array,  known  as  the  virtual 
core,  is  the  optimum  configuration.  Within  this  virtual  core 
is  found  25%  of  the  total  collection  area  of  the  LOFAR 
array.  Approximate  peak  luminosities  taken  from  Figure  8 
are  4  x  103  Jy  at  15  MHz  and  5  x  102  Jy  at  1 10  MHz.  For  a 
shock  front  with  an  angular  size  of  «30',  this  gives  ~  1  mJy 
arcscc  2  and  ~  1 00  /Jy  arcscc  2,  respectively.  At  these  inten¬ 
sities,  LOFAR  will  be  able  to  detect  any  cluster  merger 
shocks  out  to  a  distance  of  ~2000  Mpc  in  the  1 5  MHz  fre¬ 
quency  band  and  ~700  Mpc  in  the  110  MHz  frequency 
band.  These  calculations  neglect  any  absorption  due  to  the 
interstellar  medium  within  the  Milky  Way. 

We  use  results  derived  for  our  standard  system  to  esti¬ 
mate  the  detectability  of  cluster  mergers  with  EGRET  and 
GLAST?  Generalizing  the  approach  of  Dermer  & 
Schlickeiser  (2002)  to  extended  sources,  the  number  of 
background  counts  Aj,g  from  extragalactic  diffuse  7-ray 
background  photons  with  energies  greater  than  E1  detected 
within  solid  angle  AH  of  a  high-latitude  source  is  given  by 


„  (AO)A0(At)KBEW0  /  E,  \ 

aB  ~  a o~l  \E100J 


where  £joo  =  100  MeV,  KB  =  (1.72  ±  0.08)  x  10-7  photons 
cm-2  s-1  MeV-1,  and  aB  =  (2.10  ±  0.03)  (Sreekumar  et  al. 
1998).  The  on-axis  effective  area  of  the  telescope  is  approxi¬ 
mated  by  the  expression  Ao(E1/Ei00)a° ,  and  A?  is  the 
exposure-corrected  time  on-source. 

The  number  of  source  counts  produced  by  photons  with 
energies  greater  than  £/  and  detected  within  solid  angle  AH 
centered  on  the  direction  to  a  source  is  given  by 


r°°  /  e  V 

Ns(>  Ey)  £-  (A/)  A o  /  dE <j>s{E)  — 
Je ■  V^ioo/ 


1  -  exp 


-AH 


7 t02AE)  +  7 t62s 


(33) 


2  See  http:/ /www.lofar.org. 

3  See  http:/ /glast. gsfc.nasa.gov. 


Here  9d(E)  is  the  energy-dependent  angular  resolution 
(point  spread  function)  of  the  telescope,  9S  is  the  angular 
radius  of  the  source,  assumed  to  be  energy-independent, 
and  (f)s(E)  is  the  source  flux,  which  we  write  as 
HE)  =  (cj)o/Eloo)(E/Eloo)~a,‘  (photons  cm-2  s-1  MeV-1), 
with  E  in  units  of  MeV.  Writing  the  vL„  spectral  power  as 
vLv  =  L0{E'/E'm)a\  where "  E'={\+z)E,  E'm  = 
(l+c)£joo,  and  av  =  2  —  as  is  the  uLv  source  spectral 
index,  we  have 


HE) 


Lq 

47rd2(1.6  x  10-6  ergs  MeV^1)^ 


5.2  x  10-10L43 
“too 


a„— 2 


(34) 


The  luminosity  distance  dL  =  1  00/fioo  Mpc,  and  the  bright¬ 
est  luminosity  reached  at  100  MeV  is  L0  =  1043  L43  ergs  s-1, 
noting  from  the  calculations  in  Figures  6  and  7  that  Z43  «  1 . 

The  maximum  sensitivity  occurs  when  AH  = 
7 t6^{E)  «  7 rd2,  giving 


Ns{>  £7)  9*  3.3  x 


10-8  Ag  (cm2)  Af(s )L43 
(1  CIq  cxv ) 


Writing  AH  =  ird2  =  9.6  x  lO-402(°),  where  9(°)  is  the  tele¬ 
scope’s  acceptance  angle  in  degrees,  the  background  count 
rate 


Abg(>  Ej)  ^  1.65  x  10-8 


91(’)Aq  (cm2)  Ar(s) 
1.1  —  a0 


(36) 


The  sensitivity  of  detection  at  significance  n„  is  therefore 
given  by 

Ns 

n° 

-  2.6  x  10-4^43  n/Aq  (cm2)  Af(s)  / 

0(°)(1  -  «o  -  otv)  \E\qq) 

(37) 

High-confidence  identification  of  a  source  requires  that 
Ns^>  1  and  n„  >  5. 

For  the  EGRET  telescope  on  the  Compton  Gamma  Ray 
Observatory,  A0  =  1200  cm2,  a0  =  0,  and  a  nominal  observ¬ 
ing  time  is  At  =  106?6  s,  with  7,  =  1.A  source  with  a  radius 
of  0.5  Mpc  will  have  an  angular  radius  9S  ~  0/3/dioo-  The 
point-spread  function  of  EGRET  is  8d{E)  =  5° {E / EWo)x , 
which  is  larger  than  the  source  angular  size  except  for  the 
nearest  sources  or  when  observing  at  the  highest  energies, 
in  which  case  the  number  of  sources  photons  is  small. 
Substituting  into  equations  (36)  and  (37)  gives 

N,{>  Ey)  ~  404L43(£7/£ioo)^_1/[(1  -  av)d2m\  , 

n<7  -  L9v^L43(£7/£ioo)q,'+0'05/[(1  -  a„)4> 0]  ■ 

Unless  L43  ^  1 ,  we  conclude  that  it  is  unlikely  that  EGRET 
could  have  detected  nonthermal  emission  from  merging 
clusters  of  galaxies. 
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This  is  in  disagreement  with  the  estimate  of  Totani  & 
Kitayama  (2000),  who  argue  that  a  large  fraction  of  the  iso¬ 
tropic  unidentified  EGRET  sources  could  be  due  to  the 
radiation  from  merging  clusters  of  galaxies.  Totani  & 
Kitayama  (2000)  make  very  optimistic  assumptions  about 
the  hardness  of  the  spectral  index  of  injected  nonthermal 
electrons  and  in  fact  claim  that  the  merger  shocks  have  such 
large  Mach  numbers  that  the  injection  index  is  nearly  equal 
to  2.0.  As  shown  in  Figure  10,  such  hard  spectral  indexes  are 
possible  only  in  the  most  centrally  peaked  clusters,  but  the 
average  cluster  merger  shock  as  inferred  from  observations 
of  galaxy  clusters  has  a  minimum  injection  index  between 
2.2  and  2.3.  Because  the  energies  of  electrons  that  Compton 
scatter  the  CMB  radiation  to  greater  than  100  MeV  energies 
exceed  ~200  GeV,  the  steeper  injection  spectrum  reduces 
the  available  power  in  these  electrons  by  1-2  orders  of  mag¬ 
nitude,  thus  accounting  for  the  difference  between  our 
expectations  and  the  claims  of  Totani  &  Kitayama  (2000), 
Kawasaki  &  Totani  (2002),  and  Colafrancesco  (2001). 
Although  it  is  possible  that  a  few  of  the  unidentified 
EGRET  sources  are  associated  with  merging  clusters  of  gal¬ 
axies,  these  would  involve  the  less  frequent  events  involving 
collisions  of  clusters  with  masses  near  1015  M0  and  those 
with  hard  spectral  indexes  or  dark  matter  density  profiles 
with  strong  central  peaks. 

We  have  not  taken  into  account  effects  of  a  merger  tree  on 
the  evolution  of  the  nonthermal  particle  distribution 
(Gabici  &  Blasi  2003).  The  protons  that  radiate  in  the 
EGRET  and  GLAST ranges  will  accumulate  and  could  con¬ 
tribute  factors  of  2  4  enhancements  because  of  successive 
mergers,  but  not  enough  to  significantly  alter  our  conclu¬ 
sions.  Our  predictions  are  in  accord  with  a  recent  analysis  of 
EGRET  data  (Reimer  et  al.  2003;  see,  however,  Scharf  & 
Mukherjee  2002)  that  fails  to  detect  7- ray  emission  from 
galaxy  clusters  that  emit  the  most  luminous  thermal  X-ray 
emission. 

For  the  Large  Area  Detector  on  GLAST ,  A0  =  6200  cm2, 
a0  =  0.16,  and  the  on-source  observing  time  for  a  scanning 
mode  lasting  fyr  years  is  A t  =  0.2  x  3.16  x  107ryr  s,  where 
the  factor  0.2  is  the  exposure  correction.  The  point-spread 
function  for  GLAST  is  9d  =  3?5(£,/£'ioo)  ^3,  which  still 
generally  dominates  the  source  angular  size  except  at  the 
highest  energies  and  for  the  closest  sources.  From  this  we 
obtain 

Ns{>  EJ  ~  m0tyrL43(E,/Em)a"-O  S4/[(0M  -  av)d2m]  , 

na  ~  l4L43^(E,/Emr+0-3/[(OM  -  av)d2m]  . 

These  estimates  indicate  that  clusters  of  galaxies  within 
«200  Mpc  will  be  detected  with  GLAST ,  but  we  predict  far 
fewer  clusters  of  galaxies  than  predicted  by  Totani  & 
Kitayama  (2000)  for  the  reasons  discussed  above.  Note  that 
the  detection  significance  is  only  weakly  energy-dependent, 
so  that  the  criterion  for  detection  becomes  the  requirement 
of  detecting  a  few  photons,  which  is  more  easily  satisfied  at 
lower  energies. 

Our  results  also  have  relevance  to  suggestions  that  clus¬ 
ters  of  galaxies  contribute  to  the  diffuse  7-ray  background. 
Loeb  &  Waxman  (2000)  have  argued  that  Compton- 
scattered  CMB  radiation  from  relativistic  electrons  acceler¬ 
ated  by  shocks  formed  during  structure  formation  can  make 
a  significant  contribution  to  the  diffuse  extragalactic  7-ray 
background,  which  has  a  featureless  power-law  spectrum 


between  a  few  MeV  and  ~  100  GeV  (Sreekumar  et  al.  1998). 
If  this  emission  is  provided  by  cluster  merger  shocks,  then 
our  results  indicate  that  this  contribution  cannot  be  signifi¬ 
cant  for  two  reasons.  First,  the  density  distributions  inferred 
from  X-ray  observations  of  thermal  bremsstrahlung  in 
Abell  clusters  imply  an  average  minimum  injection  spectral 
index  that  would  produce  Compton-scattered  CMB  emis¬ 
sion  that  is  softer  than  the  spectrum  of  the  diffuse  7-ray 
background.  Although  those  few  clusters  with  the  hardest 
indexes  would  contribute  the  most  emission  above  100 
MeV,  the  superposition  of  the  emission  from  the  softer 
sources  would  produce  a  concave  spectrum  at  lower  ener¬ 
gies,  which  is  not  observed.  Second,  an  associated  hadronic 
signature  would  be  observed  even  in  the  case  where  protons 
are  accelerated  with  the  same  efficiency  as  the  electrons  at 
all  times  except  near  tcon  (see  Figs.  6  and  7).  The  hadronic 
feature  would,  however,  be  seen  in  the  time-averaged  7-ray 
spectrum,  so  that  calculations  omitting  hadronic  emissions 
are  incomplete  (Keshet  et  al.  2003).  If,  as  argued  in  §  3,  the 
hadronic  acceleration  efficiency  exceeds  that  of  electrons, 
then  a  7r°-decay  signature  should  be  apparent  in  the  back¬ 
ground  radiation  (see  Fig.  9)  unless  accretion  shocks  make 
the  dominant  source  of  high-energy  7-ray  cluster  emission. 
In  this  case,  most  of  the  7-ray  emission  should  be  located  on 
the  periphery  of  the  cluster,  and  the  7r0  7-rays  would  be 
reduced  dramatically  because  of  the  decreased  lower  gas 
density  in  the  proximity  of  the  protons  (see  Miniati  2003  for 
a  spatially  resolved  model  including  both  proton  and 
electron  emissions). 

Should  it  be  the  case  that  cluster  merger  shocks  do  pro¬ 
duce  a  large  fraction  of  the  diffuse  extragalactic  7-ray  back¬ 
ground,  then  both  low  hadronic  efficiencies  and  dark  matter 
halos  with  strong  central  density  peaks  in  the  majority  of 
rich  clusters  would  seem  to  be  required  (see  Fig.  10).  Under 
the  assumptions  of  isothermality  assumed  by  the  violent 
relaxation  of  dark  matter  halos  that  collapse  from  uniform 
expansion  of  the  Hubble  flow,  dark  matter  halos  should  not 
have  central  density  cusps.  Furthermore,  X-ray  temperature 
and  luminosity  profiles  are  inconsistent  with  strongly 
peaked  central  dark  matter  density  profiles  (Makino, 
Sasaki,  &  Suto  1998;  Wu  &  Chiueh  2001)  despite  the  predic¬ 
tions  of  .V-body  simulations  (Navarro,  Frenk,  &  White 
1997).  Clusters  of  galaxies  could  still  make  a  small  contribu¬ 
tion  of  order  ~1%— 10%  to  the  extragalactic  diffuse  7-ray 
background  radiation  (Colafrancesco  &  Blasi  1998;  Miniati 
2002). 

Our  simulations  show  that  protons  may  reach  ~  1 0 1 9  eV 
for  standard  parameters  but  are  not  accelerated  to  higher 
energies,  in  agreement  with  the  gyroradius  estimates  of 
Norman,  Melrose,  &  Achterberg  (1995)  for  maximum 
cosmic-ray  energies.  The  maximum  energies  of  protons  are 
primarily  constrained  by  the  available  acceleration  times  at 
early  times.  At  later  times,  the  proton  Larmor  radii  may 
exceed  the  size  scale  of  the  shock,  so  that  the  protons  then 
diffuse  into  the  surrounding  galaxy  cluster  of  galaxies  and 
undergo  no  further  acceleration.  Competition  of  the  accel¬ 
eration  and  pj  energy-loss  rates  plays  a  small  role  in  limiting 
the  highest  energy  of  protons  for  our  standard  parameters. 

For  a  diffusion  coefficient  corresponding  to  a 
Kolmogorov  spectrum  of  turbulence,  >1016  eV  protons 
can  diffuse  from  the  cluster  on  a  Hubble  time  (Berezinsky, 
Blasi,  &  Ptuskin  1997).  Assuming  Bohm  diffusion,  however, 
Berezinsky  et  al.  (1997)  and  Volk,  Aharonian,  & 
Breitschwerdt  (1996)  showed  that  clusters  of  galaxies 
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provide  a  storage  room  for  cosmic  rays  with  energies  <  1018 
eV  on  cosmological  timescales.  Our  calculations,  which 
assume  Bohm  diffusion,  confirm  that  MO18  eV  protons 
remain  within  the  cluster  for  times  much  greater  than  the 
lifetime  of  the  shock  front  (see  Fig.  4).  For  energies  >  1019 
eV,  the  loss  of  protons  from  diffusion  through  the  ICM  will 
become  significant  over  the  lifetime  of  the  shock  front  and 
will  become  a  restriction  to  the  maximum  energies.  The 
maximum  energy  imposed  by  the  Larmor  radius  size-scale 
criterion  is  always  less  than  the  maximum  energy  deter¬ 
mined  by  the  diffusion  loss  timescale  for  our  models.  For  a 
given  magnetic  field,  the  Larmor  radius  of  nonthermal  par¬ 
ticles  is  weakly  dependent  on  the  cluster  mass  (cxM1/3),  and 
despite  the  possible  range  of  masses  for  M\ ,  the  maximum 
particle  energy  Emax  ex  M 1 7 3  and  varies  little  over  the  mass 
range  of  M\.  Kang  et  al.  (1997)  argue  that  more  rapid  par¬ 
ticle  acceleration  can  take  place  in  perpendicular  shocks, 
but  this  requires  preferential  orientations  of  the  magnetic 
field  with  respect  to  the  shock  direction. 

5.  SUMMARY  AND  CONCLUSIONS 

We  have  presented  the  results  of  a  computer  code 
designed  to  calculate  the  nonthermal  particle  distributions 
of  a  shock  front  formed  in  the  merger  event  of  two  clusters 
of  galaxies.  We  have  calculated  nonthermal  particle  energy 
spectra  for  primary  electrons  and  protons  and  secondary 
electrons.  Photon  spectra  were  calculated  for  bremsstrah- 
lung,  Compton,  and  synchrotron  processes  as  well  as  7r°- 
decay  7-radiation  from  p-p  collisions.  Our  results  apply  to 
shocks  that  form  at  the  interaction  boundary  between  two 
merging  clusters  of  galaxies  and  not  to  cluster  accretion 
shocks  that  form  at  the  outer  regions  of  the  cluster. 
Evidence  for  nonthermal  particle  acceleration  in  merger 
shocks,  which  penetrate  into  the  inner  region  of  the  accret¬ 
ing  cluster,  is  provided  by  radio  halos  and  relics  and 
nonthermal  X-ray  emission  from  clusters  of  galaxies. 

The  thermal  X-ray  bremsstrahlung  was  modeled  from 
observations  of  luminosities,  temperatures,  and  masses  of 
clusters  of  galaxies.  We  modeled  the  nonthermal  emission 
under  the  assumption  that  electrons  and  protons  are  accel¬ 
erated  with  a  5%  efficiency  of  the  available  gravitational 
energy  that  is  dissipated  during  the  course  of  gravitational 
interactions  between  two  merging  clusters  of  galaxies. 
Particle  acceleration  at  shocks  formed  between  merging 
clusters  of  galaxies  produces  a  population  of  nonthermal 
electrons  that  Compton  scatters  CMB  photons.  This 
process  naturally  produces  a  power-law  distribution  of 
photons  in  the  40-80  keV  energy  range  with  luminosities 
consistent  with  HXR  observations  of  the  galaxy  clusters 


A1656,  A2256,  A3667,  and  possibly  A2199.  In  addition,  the 
EUV  emission  observed  in  A 1656  is  also  consistent  with  a 
nonthermal  Compton-scattered  CMB  origin.  Cluster 
magnetic  field  strengths  ~0. 1  //G  are  in  accord  with  this 
interpretation  of  the  observed  F1XR  excess  and  EUV 
emission. 

In  contrast  to  the  results  of  Totani  &  Kitayama  (2000), 
Kawasaki  &  Totani  (2002),  and  Colafrancesco  (2001),  we 
have  argued  that  it  is  unlikely  that  more  than  a  few  of  the 
isotropic  unidentified  EGRET  sources  can  be  attributed  to 
radiation  from  nonthermal  particles  produced  by  cluster 
merger  shocks.  Previous  studies  have  assumed  hard  non- 
thermal  spectra  that  give  brighter  7-ray  emission  than 
implied  by  our  simulation  results.  Such  hard  spectra  can  be 
obtained  only  in  the  infrequent  merger  events  involving  two 
very  massive  clusters,  or  between  clusters  where  the  dark 
matter  density  profiles  are  centrally  peaked. 

The  diffuse  extragalactic  7-ray  background  is  a  feature¬ 
less  power  law  with  photon  index  of  2.10  ±  0.03.  The  domi¬ 
nant  nonthermal  7-radiation  components  in  cluster  merger 
shocks  include  Compton-scattered  CMB  radiation  and 
bremsstrahlung  from  nonthermal  electrons  and  7r°-decay 
emission  from  nuclear  interactions  involving  nonthermal 
protons.  The  7r°-decay  signature  will  be  present  if  the 
hadronic  acceleration  efficiency  exceeds  the  electron  acceler¬ 
ation  efficiency,  as  expected  in  diffusive  shock  acceleration 
theory.  Furthermore,  spectra  calculated  using  parameters 
obtained  from  ROSAT  observations  of  45  Abell  clusters 
(Wu  &  Xue  2000)  produce  spectra  with  slopes  «2. 2-2.4. 
Unless  dark  matter  density  profiles  are  centrally  peaked  and 
hadronic  acceleration  efficiency  is  low,  our  results  imply 
that  nonthermal  emission  from  merging  clusters  of  galaxies 
can  make  only  a  minor  contribution  to  the  diffuse  extraga¬ 
lactic  7-ray  background.  Using  standard  diffusive  shock 
acceleration  with  mean  ICM  magnetic  fields  <  1  /iG,  our 
results  also  indicate  that  the  merger  shocks  in  clusters  of  gal¬ 
axies  do  not  accelerate  >1019  eV  cosmic  rays,  although 
additional  effects,  such  as  shock  obliquity  and  the  presence 
of  preexisting  particle  populations,  could  permit  higher 
energy  acceleration  by  these  shocks. 
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APPENDIX  A 

PARTICLE  ENERGY-LOSS  RATES 
Al.  ELECTRON  LOSS  RATES 

Nonthermal  electrons  lose  energy  through  Coulomb  or  ionization  interactions  with  background  thermal  electrons.  Energy 
is  lost  through  this  process  at  the  rate 


K'.coul  (  Re :  0  —  V\  |eO(  M 
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where  we  have  made  use  of  the  following  definitions: 
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The  quantity  «icm  is  the  proton  density  of  the  ICM,  8ec  is  the  electron  velocity,  Te(t)  is  the  electron  temperature  given  in 
kelvins,  and  k  is  Boltzmann’s  constant.  Nonthermal  electrons  diffuse  in  energy  space  according  to  the  diffusion  coefficient  for 
Coulomb  scattering 
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(Miller,  Larosa,  &  Moore  1996). 

For  a  fully  ionized  ICM,  the  bremsstrahlung  energy-loss  rate  (Blumenthal  &  Gould  1970)  is  given  by  the  expression 

/  8c( 


-^e,brem(Xe}  ?)  —  044  7e  T  0.36)(Ae  WleC  )  , 


(A6) 


where  r/[le  =  1.3  is  an  enhancement  factor  for  Helium  that  results  from  the  Z(Z  +  1)  dependence  of  the  loss  rate.  The  synchro¬ 
tron  energy-loss  rate  for  electrons  averaged  over  pitch  angle  in  a  magnetic  field  of  strength  B ,  Lorentz  factor  je,  and  velocity 
8e  is  (Rybicki  &  Lightman  1979) 
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where  oj  is  the  Thomson  cross  section. 

Electrons  lose  energy  by  Compton-scattering  ambient  photons  in  the  ICM  environment,  as  described  by  equations  (13), 
(16),  and  (17)  for  the  cluster  X-ray,  stellar,  and  CMB  radiation  fields,  respectively.  The  spectral  distribution  of  each  radiation 
field  is  approximated  by  a  blackbody  spectrum  of  temperature  Tx  (see  eq.  [9]),  4000  K,  and  2.7(1  +2)  K  at  epoch  z.  The 
Compton  energy-loss  rate  is  given  by 

X^xomp  (Ee ,  t)  —  —  COjye  Be  ^  (A8) 

i 

(Skibo  1993).  The  quantity  OKN(7e,  i?,-)  =  oc^d,-) l[oj{\  +  7ed;)]  is  a  Klein-Nishina  correction  factor,  where  =  kTj/mgC1  is 
the  dimensionless  temperature  and  crc(e7)  is  the  Compton  scattering  cross  section  for  photons  with  energy  e7  =  hv/mec 2  in  the 
electron  rest  frame. 

Electrons  can  also  lose  energy  through  the  double  (or  “  radiative  ”)  Compton  process  in  which  a  second  photon  is  created 
during  an  electron-photon  interaction.  Estimates  by  Ram  &  Wang  (1971)  and  Gould  (1975)  show  that  the  ratio  of  energy-loss 
rates  through  double  Compton  scattering  and  ordinary  Compton  scattering  is  when  the  typical  photon  energy  in  the 
electron  rest  frame  is  ~50  MeV.  This  ratio  only  reaches  unity  when  the  rest-frame  photon  energy  is  ~5  TeV.  Synchrotron 
energy  losses  completely  dominate  at  these  energies,  so  the  energy-loss  rate  through  double  Compton  scattering  can  be  safely 
neglected.  A  detailed  treatment  of  electron  energy  losses  in  galaxy  clusters,  although  at  energies  where  double  Compton  is  not 
important,  is  given  by  Sarazin  (1999). 


A2.  PROTON  LOSS  RATES 
The  Coulomb  energy-loss  rate  for  protons  is  given  by 
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In  addition  to  energy  losses,  the  protons  diffuse  in  energy  space  according  to  the  diffusion  coefficient 
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(Huba  1994). 
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Protons  undergo  proton-proton  collisions  to  produce  pions  in  which  a  proton  loses  of  its  energy  per  collision.  After  three 
collisions,  the  proton  is  assumed  to  be  lost  from  the  system.  The  timescale  for  this  loss  is 

Tpion(£p’ 0  =  c/VhcMK  +  O.lane)  (AU) 

(Sturner  et  al.  1997),  where  oh  and  one  are  the  inelastic  cross  sections  for  pion-producing  collisions  with  protons  and  helium, 
respectively.  Measured  values  for  these  cross  sections  are  cth  =  28  mbarns  and  <7He  =  100  mbarns,  well  above  threshold 
(Meyer  1972). 

High-energy  protons  interact  with  the  CMB  photons  because  of  photomeson  production.  It  requires  about  five  interactions 
for  protons  to  lose  a  substantial  fraction  of  their  energy  because  of  interactions  with  photons  with  proton  rest-frame  energies 
in  the  range  200-500  MeV,  and  about  two  interactions  for  photons  whose  energies  in  the  rest  frame  of  the  protons  are  greater 
than  500  MeV.  The  timescale  for  the  catastrophic  loss  is  given  by 

c  f°° 

'T~p~/—>7r(<Ep,  i)  —  —  J  derive)  J  dp  (1  p^kp^^c  )(j^_,7r(e  )  ,  (A12) 

where  apy-,v(d)  is  the  cross  section  of  the  proton-photon  interaction  and  d  =  'ype(  1  -  p)  is  the  photon  energy  in  the  rest  frame 
of  the  proton.  Here  p  =  cos  0,  and  9  is  the  angle  between  the  directions  of  the  interacting  photon  and  proton.  We  let 
(jp~,^T\d)  =  380  and  120  //barns,  and  the  inelasticity  coefficient  kpl  =  0.2  and  0.6,  for  200  MeV  <  d  <  500  MeV  and 
d  >  500  MeV,  respectively  (Atoyan&  Dermer  2001).  The  number  density  of  CMB  photons  at  redshift  zis  denoted  by  n7(e). 

Protons  diffuse  through  the  ICM  at  a  rate  that  depends  on  the  diffusion  coefficient.  This  diffusion  process  is  a  random  walk 
that  allows  particles  to  escape  from  the  system.  Protons  that  escape  from  the  cluster  merger  are  treated  as  a  catastrophic  loss 
with  timescale  (Volk  etal.  1996) 
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where  we  employ  the  Bohm  diffusion  coefficient  «: Boh m  ('£),)  =  Act'l/3  for  particle  diffusion,  and  r l  is  the  particle  Larmor 
radius. 


A3.  PION  PRODUCTION  FROM  p-p  COLLISIONS 

Nonthermal  protons  collide  with  ambient  gas  particles  to  produce  secondaries  through  several  channels.  The  five  dominant 
channels  are  (1)  p  +p  — *  7r°  +  X,  (2)  p  +p  — >  n+  +  X,  (3)  p  +p  — >  +  X,  (4)  p  +p  — >  K+  +  X,  and  (5)  p  +p  — >  K~  +  X. 

Here  X refers  to  any  by-products  produced  in  the  reaction  other  than  the  particle  indicated.  It  is  also  understood  that  channel 
2  does  not  include  the  contribution  resulting  from  the  production  of  deuterium  in  the  reaction  p  +  p  — >  7r+  +  d,  which  can  be 
treated  separately  but  is  a  small  contribution  to  the  total  secondary  production.  The  secondary  electron  production  is 
dominated  by  the  first  two  pion  channels,  and  these  are  the  only  two  channels  considered  in  our  calculations.  Our  calculation 
technique  follows  the  method  described  by  Dermer  (1986a,  1986b)  and  Moskalenko  &  Strong  (1998);  see  also  Blasi  & 
Colafrancesco  (1999).  This  model  combines  the  isobaric  model  (Stecker  1970)  at  energies  less  than  3  GeV  and  the  scaling 
model  (Badhwar,  Golden,  &  Stephens  1977)  for  energies  greater  than  7  GeV.  We  use  a  linear  combination  to  join  the  two 
models  in  the  transition  region  between  3  and  7  GeV. 

For  proton  energies  less  than  3  GeV,  the  production  of  pions  produced  by  proton-proton  collisions  is  mediated  by  the 
excitation  of  a  A3/2  isobar.  Assuming  the  outgoing  A3/2  isobar  of  mass  m,\  is  collinear  with  the  initial  direction  of  the  colliding 
protons  in  the  center-of-momentum  system  (CM)  with  isotropically  distributed  decay  products  in  the  laboratory  system  (LS), 
the  pion  distribution  in  the  LS  of  the  isobaric  model  is  given  by 


f{Kp,Kp,mA)  = 


1 


4m^7;/7(r7+/3+ 


H\id,  H+H,  H+(+)]  +  id  *r  (-),  *T  (+)] 

Ta^a 


(A  14) 


where  the  Heaviside  function  H[x;  a ,  b]  is  defined  to  be  =1  if  a  <  x  <  b,  otherwise  =0.  The  function  =  ^7^(1 

and  77  is  the  pion  LS  kinetic  energy.  The  forward  (+)  and  backward  (— )  moving  isobar  Lorentz  factors  are 
7a  =  7c7a(1  ^  PcPa)’  where  7 c=yfs/2mp  is  the  Lorentz  factor  of  the  CM  with  respect  to  the  LS,  and 
7^  =  (a  +  m\  +  mi)  j  2d  dm  A  is  the  Lorentz  factor  of  the  isobar  in  the  CM.  The  pion  Lorentz  factor  in  the  rest  frame  of  the  A 
isobar  is  7^  =  (m\  +  m~  —  mj) /2mAm~. 

The  pion  spectrum  resulting  from  the  p-p  collisions  is  calculated  by  evaluating  the  following  integral  over  the  isobar  mass 
spectrum: 


dN{Kn,Kp) 

dT^ 


Wr(Kp ) 


yfs—mp 


dmABip  (mA)f(K7[:  Kp.  mA)  . 


mp+mn 


(A15) 
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The  term  B.r(m^)  is  the  normalized  isobar  mass  spectrum  given  by  the  Breit-Wigner  distribution 

wr(Kp)T 


Bir{mK)  = 

with  average  isobar  mass  (wa),  and  normalization  factor 


0«a  -  K))-  +  r2 


wr{Kp)  = 


arctan 


( V~s -mP-  (mA)\  _  arctan  / mp  +  mw  -  (mA) 


(A  16) 


(A17) 


Our  calculations  only  consider  the  A3/2(1232)  isobar,  which  has  a  resonance  width  T  =  \  x  115  MeV  (Particle  Data  Group 
1984). 

For  protons  whose  energies  are  greater  than  7  GeV,  we  adopt  the  scaling  model  of  Stephens  &  Badhwar  (1981).  The 
Lorentz  invariant  cross  section  given  by  Badhwar  et  al.  (1977)  and  Stephens  &  Badhwar  (1981)  for  neutral  and  charged  pion 
production  is  given  by 


""  d3pn 

where  we  have  made  use  of  the  following  definitions: 


d  3  Q- 

Ek— —  =  y/^n(Ep){  1  -  exp 


-@P±  \ 
1  +  4  mils) 


9AEP) 


=  (i+h£)' 


s  ) 


/  4  m2\ 

^AEP)  =  (l  +  23Ep26)(  1--^ 


J  = 


C\  -  C2p  +  C2p2± 
\/l  +  4m2p/s 


Cr=  \  jx*  +  (-)Oi  +  ml)  , 


r*  - 

An  — 


2m^77r7fS,1/2(/37r  cos  9  -  (3C ) 
\J  (s  —  m2  —  mx)~—4m2m2x 


(A18) 

(A19) 

(A20) 

(A21) 

(A22) 

(A23) 


The  constants  -y/,  C  1,2,3,  and  Sk  are  given  by  Badhwar  et  al.  (1977)  and  Stephens  &  Badhwar  (1981),  and  0  is  the  pion  LS 
polar  angle.  Adopted  values  for  the  constants  s# ,  J?,  CT23,  and  M  are  given  in  Table  1.  The  quantities  mx  depends  on  the 
channel  under  consideration  and  are  (1)  mx  =  2 mp,  (2)  mx  =  mp  +  mn,  (3)  mx  =  2m p  +  m, r,  (4)  mx  =  mp  +  m„,  and  (5) 
mx  =  2m  f,  +  mK. 

To  calculate  the  LS  energy  distribution  of  pions,  we  integrate  the  following  over  the  pion  LS  polar  angle: 


QiriEn,  Ep) 


rl 

(r]an(Ep))  3 Cos 9max 


d  cos  9 


where  —  1  <  cos  0max  <  1  and  is  defined  as 


COS  Hmax 


1  f  s-  m2xm2\ 

PcIcPTr  v'  *  2Vs  )  ' 


(A24) 


(A25) 


The  quantity  (r]an(Ep)}  is  the  inclusive  cross  section,  and  its  values  are  taken  from  equations  (5)  (9)  of  Dermer  (1986a). 


TABLE  1 


Adopted  Values  for  Constants  Defining  the  Invariant 
Cross  Section  for p+p  —*  X  Reactions 


Channel 


Constant  1  2  3 


si .  140  153  127 

St .  5.43  5.55  5.3 

Ci .  6.1  5.3667  7.0334 

Ci .  -3.3  -3.5  -4.5 

C3  .  0.6  0.8334  1.667 

m .  2  1  3 
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The  pions  will  decay  to  produce  muons.  These  muon  are  created  fully  polarized.  This  results  in  a  e*  decay  asymmetry.  The 
muons  created  by  the  pion  decay  are  created  in  the  rest  frame  of  the  pion  with  Lorentz  factor 
In  =  {ml  +  m2)/2mwmfl  «  1.039,  and  «  0.2714.  In  the  rest  frame  of  the  muons,  the  resulting  electron-decay  spectrum  is 
given  by  (Commins  1973;  Orth  &  Buffington  1976) 


d2N 

dEf  d  cos  9* 


2e2{e  —  2e) 


1  +£(  3  _  ^  )  cos 6* 


(A26) 


where  E *  is  the  electron  energy  in  the  muon  rest  frame  and  helps  to  define  the  quantity  e  =  2E*  /mfl.  The  angle  between  the 
polarization  angle  of  the  parent  muon  and  the  resulting  electron  momentum  is  9* .  The  quantity  £  =  ±  1  for  ft+  — >  e± . 

For  a  given  electron  energy  Ee,  the  normalized  electron-decay  energy  distribution  in  the  frame  where  the  pions  are  isotropic 
is  (Dermer  1986a) 


dN(Ee;y„ ) 
dEe 


'  <t>{S-)  -  <£(«+)  : 

-  ip(u-)  +  (/){oj-)  -  <£(?+)  , 
ip{u+)  -  1p{u)-)  +  (j){u-)  -  ^(c+)  , 
#?-)  -  ip{<;+)  , 

V>(w+)  -  V'(g)  , 

-0, 


c+  <  <T-  <  0J-  , 

C+  <  u>-  <  C-  <  , 

S+  <  u>-  <  UJ+  <  C-  , 
u>-  <  s+  <  C-  <  W+  , 
<  C+  <  W+  <  C-  , 
W+  <  c-  , 


where  we  have  introduced  the  following  definitions: 


<£t  = 


hn~mii  ( 1  4:  ) 

1  ±  Pn 


u±  = 


1p{x)  = 


</>{x)  = 


1 


8  llx2 


dEg  d cos  9 * 


(3-up2){\-p2)  4x{3  +  p2-4up2) 


(5  +  u)  In  x  — 


2  9 

6  {u  +  2u/3fl  +  3)x2  16  (u  +  3up/l  +  2)x3 


(1  +  Pgf 


3(1  +Pnf 


(A27) 


(A28) 


where  u  =  £//3/(,  and  £  =  1  for  positrons  and  —  1  for  electrons.  Figure  1 1  shows  the  energy  distributions  of  electrons  and 
positrons  formed  through  p-p  collisions  by  monoenergetic  protons  with  a  variety  of  energies  for  channels  1  and  2  (compare 
Murphy,  Dermer,  &  Ramaty  1987). 

The  resulting  electron  or  positron  spectrum  formed  by  the  decay  of  isotropically  produced  pions  is 


Qe{Ee) 


dN{Ee;yn) 

dEe 


(A29) 


p  +  p  -*7T+  -»e*  p  +  p  -*»i"  -»e_ 


Fig.  1 1 . — Energy  distributions  of  secondary  positrons  (left)  and  electrons  (right)  produced  in  p-p  nuclear  collisions  by  protons  with  energies 
Kr  =  3.16  x  10s,  5.62  x  10s,  3.16  x  109,  2.15  x  1010,  1.0  x  1012,  2.5  x  1014,  6.3  x  1016,'  1.58  x  1019  ,  4.0  x  1021,  and  1.0  x  1024  eV  ( from  bottom  to  top)  for  the 
p  +p  —>  e+  +  X  channel  and  protons  of  energies  Kp  =  1.0  x  109,  2.15  x  109,  5.0  x  109,  5.0  x  1010,  1.0  x  1012,  2.5  x  1014,  6.3  x  1016,  1.58  x  1019,  4.0  x  1021, 
and  1 .0  x  1024  eV  ( from  bottom  to  top)  for  the p  +  p  — »  e~  +  X  channel. 
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with  Q^iln)  given  by  equation  (A24)  and  =  {Ee/E^m  +  E^ax/Ee)/2  if  Ee  >  E^ax  =  mM(l  +  ^ll)ryfl/2  w  69.8  MeV; 
otherwise,  %  =  1 . 


A4.  KNOCK-ON  ELECTRONS 


Nonthermal  protons  interact  with  thermal  electrons  through  Coulomb  collisions  to  produce  a  population  of  nonthermal 


knock-on  electrons.  The  knock-on  electron  source  function  in  electrons  cm~3  s~ 1  is  given  by 


Qk{le)  =  1.75«icM(r)47r  /  dip® 


'p^p^e 


(A30) 


(  Abraham,  Brunstein,  &  Cline  1966),  where  the  factor  1.75  accounts  for  the  presence  of  heavier  nuclei  in  the  ICM  and  cosmic 
rays,  7,.  is  the  Lorentz  factor  of  the  electron,  jp  is  the  Lorentz  factor  of  the  proton,  and  dJp / d^p  is  the  differential  cosmic-ray 
proton  intensity  in  particles  cm-2  s  ' 1  sr_1.  The  differential  probability  for  the  production  of  an  electron  having  a  total  energy 
per  unit  rest-mass  energy  in  d^e  at  7,.  by  the  collision  of  a  cosmic-ray  proton  of  Lorentz  factor  in  units  of  cm2  s-1  is  given  by 


^p^e{le,lp)d'ye 


2-Krj  \  1  niefap  +  (m2  +  mj)  /  (2mpme)]  m2e  \  , 

(!  -  Ip1)  \(7e  -  l)2  mp(rie  -  1)7 p  2™hp)  7e  ’ 


(A31) 


(  Abraham  et  al.  1966),  where  re  is  the  classical  radius  of  the  electron.  The  maximum  transferable  energy  is  given  by 

=  ,  _ Ip  ~  1 _ 

7max  (me/mp)  [7 p  +  {ml  +  mj)/ (2 memp)] 

The  lower  bound  of  the  integration  in  equation  (A30)  is  determined  by  solving  the  inequality  <  7max  for  jp.  The  result  is 


(A32) 


mi 


mi 


nn 


1 1/2 


(A3  3) 


APPENDIX  B 
PHOTON  PRODUCTION 

Nonthermal  photons  are  produced  by  electron  synchrotron  radiation,  electron-electron  and  electron-nucleon 
bremsstrahlung,  and  Compton-scattered  CMB,  stellar,  and  X-ray  radiation  fields.  These  components  are  calculated  for  both 
primary  electrons  and  secondary  electrons  and  positrons.  Neutral  pion-decay  7-rays  are  also  produced  in  proton-nucleon 
interactions.  We  discuss  each  component  in  the  order  mentioned  above. 


Bl.  SYNCHROTRON  RADIATION 


We  define  the  synchrotron  radiation  emissivity  as  Qsynch  =  dNsynch/dV  dtdE T  It  is  calculated  from  the  expression 

(  47t\/3(?3  b\  r11  fKl- 

Qsynch^,  t)  =  l-Tj - 3-  /  ddsm-(6)  dKeJe{Ke,t)F{is/vc)  ,  (Bl) 

V  lTvme  c  )  J0  JK„  m  m 

where  0  is  the  electron  pitch  angle,  E1  =  hi>  is  the  photon  energy  of  frequency  v,  vc  =  0.4272R_7  sin(d)  in  units  of  Hz,  and 
Je{Ee,t)  is  the  electron  intensity  (in  units  of  electrons  cm~2  s  1  sr  1  E(:v).  The  function  F{v/vc)  is  defined  through  the 
expression 

F{v/vc)  =  )  [  dxK5/3(x)  ,  (B2) 

\vcJ  Jv/Vc 

where  K5/3  is  the  modified  Bessel  function  of  order  5/3.  The  integration  is  performed  numerically  with  approximate 
expressions  for  F{v/vc )  given  by  Ginzburg  &  Syrovatskii  (1965).  Electrons  are  assumed  to  be  isotropically  distributed  in  pitch 
angle.  In  addition,  we  include  the  effects  of  free-free  absorption  of  the  synchrotron  spectra  in  the  source  using  the  absorption 
coefficient  (Ginzburg  &  Syrovatskii  1965) 


«ff  (cm  *)  =  10  2 


n1 

"ICM 

2t3/2 


17.7  + In 


(B3) 


where  the  electron  temperature  Te  is  given  in  kelvins,  the  ICM  particle  number  density  is  in  units  of  cm  3,  and  the  frequency  v 
is  in  hertz. 
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B2.  BREMSSTRAHLUNG  RADIATION 

We  calculate  the  electron-electron  (e-e)  and  electron-nucleon  (e-n)  bremsstrahlung  emissivity  by  the  expression 

rKe, mas  /  da  \ 

2brem,e-e,e-«  (-^7  i  t)  =  4tT/7](  \|  Aj je  I  <1 K rJ  1  A, .  ?)  f  j  ,  (B4) 

where  we  have  included  the  correction  factor  for  the  presence  of  helium.  Approximate  values  for  the  helium  correction 
factor  are  A'nc  =  1.2  and  A"Ie  =  1 .4.  The  expressions  for  the  differential  cross  sections  are  adopted  from  equation  (A1 )  in  Haug 
(1975)  for  electron-electron  interactions  and  from  equation  (3BN)  in  Koch  &  Motz  (1959)  for  electron-nucleon  interactions. 


B3.  COMPTON  SCATTERING 


We  calculate  the  Compton  scattering  emissivity  using  the  full  Klein-Nishina  cross  section  (Jones  1968;  Blumenthal  &  Gould 
1970)  for  relativistic  electrons  (7^  1)  by  the  expression 


,2\2 


2comp,/ (J+ ,  /)  —  Btt/'o  (/uf.c  ) 

We  define 


dE t 


T  dEeE-2Je(Ee,t) 

e t  j£;.hresh(£,) 


2t/  ln(g)  +  (1  +  2#)(1  -  9) 


(?£e,r)2(l  -q) 


2(1+  qeej) 


(B5) 


— 


4  EeE, 
(mec2)2 


q  = 


ee^Ee-Ej) 


J*resh(7q)  =  y  |  l  +  l/l 


■2  + 


(. mec 2) 
EnE, 


(B6) 

(B7) 

(B8) 


«/(£/) 


\SUjE} 


(B9) 


The  energy  E,  is  target  photon  energy,  and  /i7hrcsh  (£))  is  the  minimum  energy  required  by  an  electron  to  scatter  a  photon  from 
E,  to  energy  E1.  The  temperature,  number  density,  and  energy  density  of  the  /th  photon  component  is  given  by  7),  iy.  and  iy, 
respectively,  from  §  2.1.2. 


B4.  PION-DECAY  7-RAYS 

We  use  the  method  developed  by  Dermer  (1986b)  to  calculate  the  7r°-emissivity  from  proton-proton  collisions.  The 
7r°-emissivity  is  calculated  by 

27  0,  t)  =  47i?zicm  /  dKpJp(Kp,  t)  da^Kp)  >  (B10) 

J  j^thresh  UAjj.1) 

where  we  have  defined  the  neutral  pion  kinetic  energy  as  Kj>,  and  the  differential  cross  section  for  7r°-production  is  calculated 
from  the  cross  section  given  by  Dermer  (1986b),  which  uses  the  isobar  model  of  Stecker  (1973)  for  proton  energies  less  than 
3  GeV  and  the  scaling  model  of  Stephens  &  Badhwar  (1981)  for  protons  of  energies  greater  than  7  GeV.  A  linear  combination 
of  the  two  models  is  assumed  between  3  and  7  GeV.  The  7-ray  emissivity  from  7r0-decay  is  calculated  from  the  expression 


/»00  QF  F 

Qp/(E1,t)=2  dE^^fdL 

JE1+my/4Ey  2o 


"Vc4 


(B 1 1 ) 


Considering  collisions  with  heavy  nuclei,  Dermer  (1986b)  has  shown  that  the  7-ray  emissivity  is  enhanced  by  the  multiplicative 
factor  1 .45.  The  total  7-ray  emissivity  therefore  becomes 


QlT' 

^~"7ru  sy 


1.452/ 


(B12) 


REFERENCES 


Abraham,  P.  B.,  Brunstein,  K.  A.,  &  Cline,  T.  L.  1966,  Phys.  Rev.,  150, 
1088 

Arnaud,  M„  &  Evrard,  A.  E.  1999,  MNRAS,  305, 631 
Atoyan,  A.,  &  Dernier,  C.  D.  2001,  Phys.  Rev.  Lett.,  87,  221102 


Atoyan,  A.  M„  &  Volk,  H.  J.  2000,  ApJ,  535, 45 

Badhwar,  G.  D.,  Golden,  R.  L.,  &  Stephens,  S.  A.  1977,  Phys.  Rev.  D,  15, 
820 

Bagchi,  J.  2003,  MNRAS,  in  press 


No.  2,  2003 


NONTHERMAL  PARTICLES  AND  RADIATION 


731 


Baring,  M.  G.,  Ellison,  D.  C.,  Reynolds,  S.  P..  Grenier,  I.  A.,  &  Goret,  P. 

1999,  ApJ,  513,311 
Bell,  A.  R.  1978,  MNRAS,  182,  147 

Berezinsky,  V.  S„  Blasi,  P.,  &  Ptuskin,  V.  S.  1997,  ApJ,  487,  529 
Berezinskii,  V.  S„  &  Grigor’eva,  S.  I.  1988,  A&A,  199,  1 
Berghofer,  T.  W„  &  Bowyer,  S.  2002,  ApJ,  565,  L17 
Berghofer,  T.  W.,  Bowyer,  S.,  &  Korpela,  E.  2000a,  ApJ,  535,  615 

- .  2000b,  ApJ,  545, 695 

Berrington,  R.  C.  2000,  Ph.D.  thesis,  Indiana  Univ. 

Berrington,  R.  C„  Lugger,  P.  M.,  &  Colin,  H.  N.  2002,  AJ,  123,  2261 

Bertschinger,  E.  1985,  ApJS,  58,  39 

Blandford,  R.  D.,  &  Ostriker,  J.  P.  1978,  ApJ,  221,  L29 

Blasi,  P.  2000,  ApJ,  532,  L9 

Blasi,  P.,  &  Colafrancesco,  S.  1999,  Astropart.  Phys.,  12, 169 
Blumenthal,  G.  R„  &  Gould.  R.  J.  1970,  Rev.  Mod.  Phys.,  42, 237 
Bonamente,  M..  Lieu,  R.,  &  Mittaz,  J.  P.  D.  2001,  ApJ,  547,  L7 
Bowyer,  S.,  &  Berghofer,  T.  W.  1998,  ApJ,  506,  502 
Bowyer,  S.,  Korpela,  E.,  &  Berghofer,  T.  2001,  ApJ,  548,  L135 
Bowyer,  S.,  Lampton,  M.,  &  Lieu,  R.  1996,  Science,  274,  1338 
Bridle,  A.  H.,  &  Fomalont,  E.  B.  1976,  A&A,  52,  107 
Brunetti,  G.,  Setti,  G.,  Feretti,  L.,  &  Giovannini,  G.  2001,  MNRAS,  320. 
365 

Bryan,  G.  L„  &  Norman,  M.  L.  1998,  ApJ,  495,  80 
Clarke,  T.  E„  Kronberg,  P.  P„  &  Bohringer,  H.  2001,  ApJ,  547,  LI  1 1 
Colafrancesco,  S.  2001.  in  AIP  Conf.  Proc.  587,  Gamma  2001:  Gamma- 
Ray  Astrophysics  2001,  ed.  S.  Ritz,  N.  Gehrels,  &  C.  R.  Shrader 
(New  York:  AIP),  427 

Colafrancesco,  S.,  &  Blasi,  P.  1998,  Astropart.  Phys.,  9, 227 

Colless,  M.  1989,  MNRAS,  237,  799 

Colless,  M„  &  Dunn,  A.  M.  1996,  ApJ,  458, 435 

Commins,  E.  1973,  Weak  Interactions  (New  York:  McGraw-Hill) 

Dernier,  C.  D.  1986a,  ApJ,  307, 47 

- .  1986b,  A&A,  157,223 

Dermer,  C.  D.,  &  Schlickeiser,  R.  2002,  ApJ,  575,  667 
Eilek,  J.  1999,  in  Proc.  Workshop  Diffuse  Thermal  and  Relativistic  Plasma 
in  Galaxy  Clusters,  ed.  FI.  Bohringer,  L.  Feretti.  P.  Schuecker,  &  V.  K. 
Kapahi  (MPE  Rep.  271;  Garching:  MPI),  71 
Ensslin,  T.  A.,  &  Biermann,  P.  L.  1998,  A&A,  330, 90 
Ensslin,  T.  A.,  Biermann,  P.  L„  Klein,  U„  &  Kohle,  S.  1998,  A&A,  332,  395 
Ensslin,  T.  A.,  &  Briiggen,  M.  2002,  MNRAS,  331.  101 1 
Ensslin,  T.  A.,  &  Gopal-Krishna.  2001,  A&A,  366,  26 
Ensslin,  T.  A.,  Simon.  P.,  Biermann,  P.  L.,  Klein,  U.,  Kohle,  S.,  Kronberg, 
P.  P„  &  Mack,  K.-H.  2001,  ApJ,  549,  L39 
Feretti,  L.,  Brunetti,  G.,  Giovannini,  G.,  Govoni,  F.,  &  Setti,  G.  2000,  in 
Constructing  the  Universe  with  Clusters  of  Galaxies,  ed.  F.  Durret  & 
D  Gerbal  (Paris:  IAP) 

Forman,  W.,  Bechtold,  J.,  Blair,  W.,  Giacconi,  R.,  Van  Speybroeck,  L.,  & 
Jones,  C.  1981,  ApJ,  243,  L133 
Fujita,  Y„  &  Sarazin,  C.  L.  2001,  ApJ,  563,  660 

Fusco-Femiano,  R.,  Dal  Fiume,  D..  Feretti,  L.,  Giovannini,  G.,  Grandi, 
P..  Matt,  G„  Molendi,  S„  &  Santangelo,  A.  1999,  ApJ,  513,  L21 
Fusco-Femiano,  R.,  Dal  Fiume,  D.,  Orlandini,  M.,  Brunetti,  G.,  Feretti, 
L„  &  Giovannini,  G.  2001,  ApJ,  552,  L97 
Fusco-Femiano,  R.,  Dal  Fiume,  D.,  Orlandini.  M.,  De  Grandi,  S., 
Molendi,  S.,  Feretti,  L.,  Grandi,  P.,  &  Giovannini,  G.  2003a,  in  ASP 
Conf.  Ser.  301,  Matter  and  Energy  in  Clusters  of  Galaxies,  ed.  S.  Bowyer 
&  C.-Y.  Hwang  (San  Francisco:  ASP),  in  press 
Fusco-Femiano,  R.,  Orlandini.  M.,  De  Grandi,  S.,  Molendi,  S.,  Feretti,  L., 
Giovannini,  G.,  Bacchi,  M.,  &  Govoni,  F.  2003b,  A&A,  398, 441 
Fusco-Femiano,  R.,  et  al.  2000,  ApJ,  534,  L7 
Gabici,  S„  &  Blasi,  P.  2003,  ApJ,  583,  695 
Geller,  M.  J.,  &  Beers,  T.  C.  1982,  PASP.  94, 421 
Ginzburg,  V.  L.,  &  Syrovatskii,  S.  I.  1965,  ARA&A,  3,  297 
Giovannini,  G.,  &  Feretti,  L.  2000,  NewA,  5,  335 
Gould,  R.  J.  1975,  ApJ,  196,  689 

Govoni,  F.,  Feretti.  L.,  Giovannini.  G.,  Bohringer,  H.,  Reiprich,  T.  H..  & 
Murgia,  M.  2001,  A&A,  376,  803 
Haug,  E.  1975,  Z.  Naturforsch.,  30a,  1099 

Huba,  J.  D.  1994,  NRL  Plasma  Formulary  (Publ.  6790-94-265; 

Washington,  DC:  NRL) 

Hwang,  C.-Y.  1997,  Science,  278,  1917 
Jones,  F.  C.  1968,  Phys.  Rev.,  167,  1159 

Kaastra,  J.  S.,  Lieu,  R.,  Mittaz,  J.  P.  D.,  Bleeker,  J.  A.  M.,  Mewe,  R., 
Colafrancesco,  S.,  &  Lockman,  F.  J.  1999,  ApJ,  519,  LI  19 
Kang,  H„  Rachen,  J.  P„  &  Biermann,  P.  L.  1997,  MNRAS,  286,  257 


Kassim,  N.  E.,  Clarke,  T.  E.,  Ensslin,  T.  A.,  Cohen,  A.  S.,  &  Neumann, 
D.  M.  2001,  ApJ,  559,  785 
Kawasaki,  W„  &  Totani,  T.  2002,  ApJ,  576,  679 

Keshet,  U.,  Waxman,  E.,  Loeb,  A.,  Springel,  V.,  &  Hernquist,  L.  2003, 
ApJ,  585,  128 
King,  I.  R.  1966,  AJ,  71,64 

Koch,  H.  W.,  &  Motz,  J.  W.  1959,  Rev.  Mod.  Phys.,  31, 920 
Large,  M.  I.,  Mathewson,  D.  S.,  &  Haslam,  C.  G.  T.  1959,  Nature,  183, 
1663 

Liang,  H..  Dogiel,  V.  A.,  &  Birkinshaw,  M.  2002,  MNRAS,  337,  567 
Lieu,  R.,  Bonamente,  M.,  &  Mittaz,  J.  P.  D.  1999a,  ApJ,  517,  L91 
Lieu,  R.,  Ip,  W.-H.,  Axford,  W.  I.,  &  Bonamente,  M.  1999b,  ApJ,  510, 
L25 

Lieu,  R.,  Mittaz,  J.  P.  D.,  Bowyer,  S.,  Breen,  J.  O.,  Lockman,  F.  J., 
Murphy,  E.  M.,  &  Hwang,  C.-Y.  1996a,  Science,  274,  1335 
Lieu,  R.,  Mittaz,  J.  P.  D„  Bowyer,  S.,  Lockman,  F.  J.,  Hwang,  C.-Y.,  & 
Schmitt,  J.  H.  M.  M.  1996b,  ApJ,  458,  L5 
Loeb,  A.,  &  Waxman,  E.  2000,  Nature,  405,  156 
Makino,  N„  Sasaki,  S„  &  Suto,  Y.  1998,  ApJ,  497,  555 
Markevitch,  M.,  &  Vikhlinin,  A.  2001,  ApJ,  563,  95 
Meyer,  J.  P.  1972,  A&AS,  7, 417 
Miniati,  F.  2002,  MNRAS,  337,  199 
- .  2003,  MNRAS,  in  press 

Miniati,  F„  Jones,  T.  W„  Kang,  H„  &  Ryu,  D.  2001,  ApJ,  562, 233 
Miniati,  F.,  Ryu,  D.,  Kang,  H.,  Jones,  T.  W.,  Cen,  R.,  &  Ostriker,  J.  P. 
2000,  ApJ,  542,  608 

Miller.  .1.  A.,  Larosa,  T.  N„  &  Moore,  R.  L.  1996,  ApJ,  461, 445 
Mittaz,  J.  P.  D.,  Lieu,  R.,  &  Lockman,  F.  J.  1998,  ApJ,  498,  L17 
Moskalenko,  I.  V.,  &  Strong,  A.  W.  1998,  ApJ,  493,  694 
Murphy,  R.  J.,  Dermer,  C.  D.,  &  Ramaty,  R.  1987,  ApJS,  63,  721 
Navarro,  J.  F„  Frenk,  C.  S„  &  White,  S.  D.  M.  1997,  ApJ,  490, 493 
Norman,  C.  A.,  Melrose,  D.  B.,  &  Achterberg,  A.  1995,  ApJ,  454,  60 
Oegerle,  W.  R„  &  Hill,  J.  M.  1994,  AJ,  107,  857 
- .  2001,  AJ,  122,2858 

Ohno,  H„  Takizawa,  M„  &  Shibata,  S.  2002,  ApJ,  577,  658 
Orth,  C.  D„  &  Buffington,  A.  1976,  ApJ,  206,  312 
Owen,  F.  N.,  &  Eilek,  J.  A.  1998,  ApJ,  493,  73 

Paolillo,  M.,  Andreon,  S.,  Longo,  G.,  Puddu,  E.,  Gal,  R.  R.,  Scaramella, 
R„  Djorgovski,  S.  G„  &  de  Carvalho,  R.  2001,  A&A,  367,  59 
Particle  Data  Group.  1984,  Rev.  Mod.  Phys.,  56,  S232 
Peebles,  P.  J.  E.  1993,  Principles  of  Physical  Cosmology  (Princeton: 

Princeton  Univ.  Press) 

Petrosian,  V.  2001,  ApJ,  557,  560 

Pierce,  M.  J.,  &  Berrington,  R.  C.  2003,  ApJ,  in  press 

Ram,  M.,  &  Wang,  P.  Y.  1971,  Phys.  Rev.  Lett.,  26, 476  (erratum  26,  1210) 

Reimer,  O.,  Pohl,  M.,  Sreekumar,  P.,  &  Mattox,  J.  R.  2003,  ApJ,  588,  155 

Rephaeli,  Y.  1979,  ApJ,  227,  364 

Rephaeli,  Y„  Gruber,  D.,  &  Blanco,  P.  1999,  ApJ,  51 1,  L21 

Reynolds,  S.  P.  1996,  ApJ,  459,  L13 

Ricker,  P.  M.  1998,  ApJ,  496,  670 

Ricker,  P.  M„  &  Sarazin,  C.  L.  2001,  ApJ,  561,  621 

Rybicki,  G.  B.,  &  Lightman,  A.  P.  1979,  Radiative  Processes  in 
Astrophysics  (New  York:  Wiley) 

Sarazin,  C.  L.  1999,  ApJ,  520,  529 

Scharf,  C.  A„  &  Mukherjee,  R.  2002,  ApJ,  580,  154 

Schindler,  S„  &  Muller,  E.  1993,  A&A,  272,  137 

Schlickeiser,  R.,  Sievers,  A.,  &  Thiemann,  H.  1987,  A&A,  182,21 

Skibo,  J.  G.  1993,  Ph.D.  thesis,  Univ.  Maryland 

Slee,  O.  B.,  Rov,  A.  L.,  Murgia,  M.,  Andernach,  H.,  &  Ehle,  M.  2001,  AJ, 
122, 1172 

Sreekumar,  P.,  et  al.  1998,  ApJ,  494,  523 
Stecker,  F.  W.  1970,  Ap&SS,  6,  377 
- .  1973,  ApJ,  185,499 

Stephens,  S.  A.,  &  Badhwar,  G.  D.  1981,  Ap&SS,  76,  213 
Sturner,  S.  J.,  Skibo,  J.  G.,  Dermer,  C.  D.,  &  Mattox,  J.  R.  1997,  ApJ,  490, 
619 

Takizawa,  M.  2000,  ApJ,  532,  183 

Totani,  T.,  &  Kitayama,  T.  2000,  ApJ,  545,  572 

Tribble,  P.  C.  1993,  MNRAS,  263,  31 

Valtchanov,  I.,  Murphy,  T.,  Pierre,  M.,  Hunstead,  R.,  &  Lemonon,  L. 
2002,  A&A,  392,  795 

Volk,  H.  J.,  Aharonian,  F.  A.,  &  Breitschwerdt,  D.  1996,  Space  Sci.  Rev., 
75,279 

Wu,  X.-P.,  &  Chiueh,  T.  2001,  ApJ,  547,  82 
Wu,  X.-P.,  &  Xue,  Y.-J.  2000,  ApJ,  542,  578 


